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FOREWORD 


Reliability  and  service  life  prediction  has  become  an  important  phase  in 
the  research  and  development  efforts  of  the  propulsion  community  and  is  vital 
to  the  logistics  and  procurement  schedule  of  missiles.  Unfortunately,  pre- 
dictive techniques  are  very  limited  and  only  crude  predictions  can  be  made. 
Several  reasons  fer  the  backward  state  of  affairs  include  the  statistical 
nature  of  the  propellant  properties,  the  variability  of  the  environment,  the 
translation  of  the  various  environmental  factors  into  engineeringly  meaning- 
ful loads,  and  the  understanding  of  the  process  of  propellant  failure. 

The  work  described  in  this  report  was  undertaken  under  an  NWC  contract 
to  Aerojet  Solid  Propulsion  Company  under  Contract  Number  N00123-72-C-2293 
as  part  of  the  Naval  Air  Systems  Command  Task  F31332301,  on  Mechanical 
Integrity  of  Solid  Propellant  Rocket  Motors,  to  the  Naval  Weapons  Center.  The 
work  described  in  this  report  also  augments  and  supplements  work  being  per- 
formed in  a fundamental  research  under  the  Naval  Ordnance  Systems  Command 
Task,  R0240201,  Physical  Chemistry  of  Filled  and  Unfilled  Elastomers.  Thic 
work  was  undertaken  in  order  to  develop  a better  understanding  of  the  failure 
process  in  solid  rocket  motors  and  to  examine  methods  by  which  the  known  non- 
linear behavior  of  propellant  in  failure  can  be  described.  This  is  important 
In  developing  useful  engineering  tools  for  the  prediction  of  failure.  Under 
this  program  a new  theory  that  is  potentially  useful  and  amendable  to  an 
engineering  prediction  was  developed.  In  addition  a number  of  nonlinear 
effects  that  have  an  effect  upon  the  predicted  stress  state,  failure  propaga- 
tion were  discovered. 

This  report  is  a facsimile  of  the  report  prepared  by  Aerojet  Solid 
Propulsion  Company. 

Released  by  Under  authority  of 

WILLIAM  S.  McEWAN,  Head  HUGH  W.  HUNTER,  Head 

Chemiatry  Division  Research  Department 

6 August  1974 
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PREFACE 

This  report  summarizes  the  experimental  efforts  conducted  by 
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and  the  theoretical  studies  of  Dr.  Paul  J.  Blatz  at  Shock  Hydrodynamics. 
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INTRODUCTION 

Non-linear  fracture  mechanics*  Is  a newly  developing  analytical 
method  for  the  treatment  of  fracture  behaviors  In  composite  solid  propel- 
lants. It  stands  in  contrast  to  the  well  known  linear-elastic  fracture 
mechanics  of  Griffith  (2)  and  the  linear-viscoelastic  fracture  mechanics 
of  Williams  (3).  These  latter  theories  assume  progressive  fracturing  of 
a single  crack,  with  small  strains  in  the  body  away  from  the  crack  tip. 

The  non-linear  theory  includes  non-linear  viscoelasticity  with  its  perma- 
nent memory  effects,  non-progressive  fracturing  of  cracks,  and  finite 
strains  away  from  the  crack  tip. 

The  study  presented  here  was  designed  to  illustrate  the  inadequacies 
of  the  linear  theories  and  the  fundamental  requirements  of  the  non-linear 
one.  An  empirically  derived  relation  forms  the  fundamental  point  of  depar- 
ture for  the  non-linear  fracture  criterion.  In  addition,  a non-linear 
strain  energy  function  was  investigated  as  the  initial  effort  in  developing 
the  required  analytical  relations. 

The  objectives  of  the  following  report  are: 

1.  To  define  the  key  characteristics  in  the  fracture  processes  of 
solid  propellants. 

2.  To  summarize  existing  mathematical  concepts  for  non-linear 
material  response  and  failure. 

3.  To  make  recommendation  for  future  studies. 


♦Suggested  as  a contrast  to  G.  R.  Irwin's  article,  "Linear  Fracture  Mechanics 
in  Relation  to  Visco-elastic  Materials"  (1). 
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SUMMARY 

The  empirical  observations  and  the  practical  relationships  derived 
from  them  represent  the  majority  of  the  report.  The  last  part  of  the 
report  covers  a newly  developed  theory  of  propellant  response. 

Empirical  Fracture  Studies 

Fracture  Initiation  and  the  trajectory  of  crack  propagation  in  solid 
propellants  are  important  to  the  prediction  of  motor  service  life  (4  and  5). 
However,  past  experience  with  solid  propellants  has  given  ample  evidence 
that  solid  propellants  do  not  follow  the  classical  relations.  Partially 
recognizing  this,  Fulbright  and  Miller  (6)  have  recommended  the  dissection 
of  motors  to  complete  their  failure  analyses. 

At  ASPC,  the  first  challenge  to  the  linear-elastic  concept  (Griffith) 
involved  a problem  of  mixed  boundaries.  Here,  the  classical  criteria  require 
a unique  relationship  between  displacements  of  the  crack  boundary  and  those 
at  the  specimen  boundary.  But,  it  is  shown  that  the  shape  of  the  crack  tip 
has  a very  strong  dependence  upon  the  past  loading  history  of  the  test  specimen. 
Thus,  there  is  not  a unique  relation  between  the  displacements  at  the  crack 
and  the  specimen  boundaries. 

As  an  added  demonstration  of  this  effect,  it  is  shown  that  notching 
a uniaxial  tensile  specimen  may  actually  increase  its  apparent  strength. 

This  effect  is  attributed  to  the  Increased  strain  rate  at  the  crack  tip. 

In  a later  section  of  this  report,  it  is  shown  that  a major  modification 
of  the  classical  criterion  can  be  used  to  fit  solid  propellant  fractures.  The 
observed  failure  relation  is: 

kr  V2  1/2 

• «.  - (1) 

This  relation  holds  for  creep-to-f allure  tests  and  states  that  the 
maximum  principal  stress  difference,  oUf)  -a«,  is  proportional  to  the  square 

1 1 /2 

root  of  the  modulus  difference,  [E(Wq)-Ee]  ' Equation  (1)  reduces  to 
Griffith's  criterion  when  the  propellant  is  glassy. 

The  free  surface  energy,  r.  Is  a constant  that  holds  over  the  entire 
range  from  glassy  to  rubbery  behaviors. 

A unique  feature  of  Equation  (1)  is  the  constant,  a,  that  is  used  to 
define  the  "effective-time"  of  the  relaxation  modulus,  E(Cf/q).  This  constant 
causes  the  entire  plot  of  the  relaxation  modulus  to  be  shifted  along  the  time 
axis. 
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All  modes  of  testing  and  past  loading  histories  can  be  made  to  fit 
Equation  (1)  through  either  linear  (7  to  10)  or  non-linear  (11  and  12) 
cumulative  damage  analyses;  whichever  applies  to  the  material  behavior. 


Thus,  Equation  (1),  combined  with  cumulative  damage  analyses,  provides 
an  empirical  basis  for  the  required  non-linear  fracture  criterion.  Note: 
this  relation  Is  not  considered  to  be  the  final  one,  but  It  will  be  employed 
until  a better  criterion  is  established. 

Verifications  of  Equation  (1)  Involved  constant  load  and  constant 
rate  failure  testing  of  specimen*  with  various  crack  depths,  C.  Both 
tests  gave  the  required  agreement,  although  the  constant  load  tests  evinced 
specimen  distortions  at  the  higher  loads  and  crack  depths. 

A number  of  factors  must  be  taken  into  account  in  developing  the 
relevant  non-linear  fracture  criterion.  These  are  listed  In  a special 
section  of  this  report. 

The  basic  mechanism  of  fracture  in  propellant  binders  seems  to  be 
initiation  in  tension,  with  propagation  occurring  under  the  peeling  action 
of  combined  tensile  and  shear  deformations. 

Crack  propagation  mechanisms  in  propellants  vary  considerably  on 
going  from  low  to  high  rates.  At  low  propagation  rates  a great  number  of 
Initiation  sites  are  generated,  with  cracks  emanating  from  each  of  them. 
However,  unlike  the  classical  process,  the  cracks  deviate  from  their 
expected  paths  of  propagation  following  curved  trajectories.  As  the  cracks 
turn,  the  rate  of  tearing  diminishes  and  the  crack  eventually  stops  Later, 
the  individual  cracks  are  joined  (rather  raggedly)  by  separate  tearing 
actions  between  them,  leading  to  an  overall  gross  crack.  A similar  process, 
but  with  special  differences,  also  occurs  at  the  propellant/liner  bond  inter- 
face. 

At  high  propagation  rates,  where  the  propellant  is  glassy,  the  classical 
fracture  processes  are  obtained. 

Thus,  the  rates  of  crack  propagation  determine  the  number  of  initiation 
sites,  the  extent  of  energy  dissipation,  the  trajectory  of  the  moving  crack 
tip,  and  the  appearance  of  the  fractured  surfaces. 

A primary  point,  of  relevance  to  Equation  (1),  is  the  observation  that 
dewetting  effects  define  the  mode  of  propellent  fractures.  It  is  observed  that 
when  at  least  one  loading  axis  undergoes  tensile  deformation  the  failure  mode 
will  be  a uniaxial  one  and  will  follow  a maximum  principal  stress  failure 
criterion. 
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Even  shear  failures  were  found  to  be  accompanied  by  dewetting  effects 
of  the  type  seen  in  uniaxial  tensile  testing.  The  dewetting  patterns  were 
observed  to  be  rotated  with  respect  to  the  specimen  axis,  as  expected. 

Loading-path  dependence,  was  demonstrated  for  three  different  loading 
histories.  One  involved  thermal -strain  paths  and  the  remaining  two  were 
variations  of  the  Bauschinger  effect  (14). 

In  the  first  case,  it  is  shown  that  different  failure  elongations 
will  be  obtained  depending  upon  the  order  of  straining  the  specimen  and 
varying  the  environmental  temperature. 

The  Bauschinger  effect  entails  a temporary  shift  of  the  stress-strain 
coordinate  axes.  In  propellants  this  "temporary  set"  is  readily  accomplished. 
All  that  needs  to  be  done  is  to  hold  the  specimen  under  a deformation  for  a 
period  of  time.  This  produces  marked  stress  relaxation  effects  that  are  seen 
as  major  shifts  In  the  coordinate  stress-strain  axes. 

The  first  tests  measured  the  effects  of  putting  a specimen  Into  5% 
compressive  strain  prior  to  crack-failure  testing,  in  notched  strip-biaxial 
tensile  specimens.  The  previous  load  history  gave  a 24%  decrease  in  the 
failure  deformation  and  a 15%  decrease  in  the  average  failure  stress. 

The  second  set  of  tests  Involved  a pre-loading  in  stress  relaxation 
under  a 10%  tensile  strain.  Then,  without  allowing  recovery,  the  specimens 
were  tested-to-failure  under  a superimposed  pressure.  The  resulting  improve- 
ment in  the  failure  envelope  showed  that  the  shift  in  the  strain  coordinate 
was  considerable,  up  to  the  full -extent  of  the  pre-strain. 

This  path  dependency  of  solid  propellants  obviously  necessitates  stress- 
strain  constitutive  relations  that:  (1)  Have  a permanent  memory;  (2)  allow  for 

propellant  dewetting;  and  (3)  include  finite  deformations.  The  theoretical 
studies  of  Blatz  are  included  in  this  report  as  one  approach  to  the  required 
response  analyses.  A separate  approach  by  Farris  (12  and  13)  Is  directed  at 
the  same  goal . 

The  Importance  of  dewetting  to  solid  propellant  failure  mechanisms  has 
been  considered  by  Lindsey  (15),  Robinson  (16)  and  Knauss  (17).  Although  we 
agree  In  principle, we  find  that  the  exact  role  of  dewetting  Is  not  e’ear. 

This  is  best  seen  from  a series  of  tests  with  several  different  HTPB  binder 
formulations.  It  was  found  that  the  tensile  strength  was  Independent  of 
filler  content;  at  least,  for  filler  fractions  varying  from  50  to  76.3  vol.  %. 

A few  data  points  Indicate  that  this  Independence  of  filler-content  may  extend 
to  the  unfilled  binder. 
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We  also  considered  binder  crosslinking  and  showed  that  propellant 
binders  form  a special  class  that  distinguishes  them  from  the  well  known 
elastomers.  In  conventional  rubbers  the  crosslinking  Is  much  greater 
than  that  used  In  solid  propellant  binders.  The  works  of  Bueche  and 
Dudek  (18,  19  and  20)  have  shown  that  failure  In  highly  crossllnked 
elastomers  exhibits  little  time-dependence,  while  the  very  lightly  cross- 
linked  propellant  binders  should  exhibit  considerable  time,  or  strain  rate, 
dependence.  Similarly,  their  data  suggest  that  propellant  failure  behaviors 
should  be  nearly  Independent  of  binder  crosslinking. 

In  the  characterization  of  solid  propellant  fracture  processes  we 
find  a number  of  difficulties.  It  is  shown,  for  example,  that  fracture 
initiation  testing  can  be  associated  with  a large  radius  of  curvature  at 
the  crack  tip.  This  becomes  great  enough  in  some  cases  to  reduce  the  local 
stress  field,  causing  failure  Initiation  to  occur  at  a point  away  from  the 
crack  tip. 

Equation  (1)  requires  a careful  evaluation  of  the  relaxation  moduli. 

But,  as  noted  by  Francis  (21),  the  relaxation  modulus  is  strongly  affected 
by  the  strain  level.  Also,  Adicoff  and  Lepie  (22)  showed  that  dewetting 
under  strain  causes  the  propellant  to  depart  from  thermorheologlcal  simplicity. 

Thus,  ws  concluded  that  the  relaxation  moduli  be  obtained  at  the  lowest 
possible  strain  levels  (l.e.  below  about  0.25%).  Here  the  moduli  values  are 
the  largest  and  the  effects  noted  by  Adicoff  and  Lepie  are  minimized. 

New  Theory  of  Propellant  Response  Behavior 

A non-linear  strain  energy  function  has  been  proposed  to  represent 
both  the  response  behavior  and  fracture  behavior  of  a CTPB  propellant  In 
various  stress  fields  and  under  superimposed  hydrostatic  pressures.  The 
strain  energy  function  includes  two  features  which  are  novel  and  useful. 

First,  a generalized  strain  measure  proposed  by  Seth  (22)  Is  used  to  repre- 
sent elastic  non-linearity.  Secondly,  to  account  for  load-unload  hysteresis, 
the  elastic  parameters  have  been  allowed  to  depend  upon  an  additional  field 
variable,  which  is  the  area  fraction  cf  dewetted  solids;  this  field  variable 
is  controlled  by  the  Griffith  power  balance.  The  theory  in  its  present  form 
contains  five  parameters. 

The  theory  has  been  checked  on  Aerojet  propellant,  on  simple  tensile 
data  at  77°F  from  0.05  min~^  to  1.25  min"'  and  under  superimposed  pressures 
from  0 pslg  to  500  psig.  The  theory  provides  excellent  correlation  with 
the  loading  data. 

In  a general  homogeneous  stress  field,  as  experienced  in  a motor,  the 
dewetting  field  will  be  non-uniform  and  the  extent  of  dewetting  will  depend 
on  the  entire  load-unload  history.  The  ultimate  goal  of  this  program  is  to 
develop  a finite  element  program,  into  which  Is  coupled  the  Griffith  criterion, 
and  which  will  predict  the  instantaneous  state  of  dewetting  of  a rocket  motor 
at  any  time  during  its  thermomechanical  history. 
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CHALLENGING  GRIFFITH 

There  seem  to  be  many  circumstances  where  the  Griffith  Failure 
Criterion  Is  not  obeyed  by  solid  propellants.  The  following  discussion 
considers  only  one  of  these  circumstances;  namely*  the  mixed  boundary 
assumption  of  the  Griffith  criterion.  A particular  example  of  the  mixed 
boundary  problem  is  shown  by  the  effects  of  strain  rate  at  the  crack  tip. 

Mixed  Boundary  Problem 

One  of  the  fundamental  features  of  Griffith's  Criterion  (2)  is  the 
definition  of  the  conditions  of  stress,  stress  Intensity*  or  energy 
that  exists  at  the  crack  tip.  An  unobtrusive  assumption  of  these  analyses 
is  that  there  is  a unique  relation  between  the  geometry  of  the  crack  tip 
and  the  loadings  imposed  on  the  boundaries  of  the  overall  specimen.  Our 
experience  with  rubbers  and  solid  propellants  ciearly  shows  this  to  be  a 
mixed  boundary  problem  where  the  loadings  on  the  crack  boundary  are  not 
uniquely  defined  by  those  on  the~extemr. 

This  problem  may  be  simply  reviewed  in  terms  of  the  Griffith  model 
shown  in  Figure  1.  Here,  a crack  of  length  2C  is  included  in  a simple 
tensile  specimen  that  is  loaded  only  along  its  upper  and  lower  boundaries. 

The  Griffith  criterion  specifically  defines  the  maximum  stress 
at  the  crack  tiD,  c,„  . in  terms  of  as  follows; 

mtSX  0 

°mx  =(°o[1  * 2 V^J)  (2) 

Where  p is  the  radius  of  the  crack  tip 

0 is  the  stress  applied  along  the  boundaries  of  the 
0 uniaxial  specimen 

Variations  in  p form  the  basis  fcr  the  non-singular  boundary  relation- 
ships mentioned  above. 

Since  p is  usually  an  experimentally  inaccessible  quantity,  the 
theoreticians  have  generally  used  classical  elasticity,  or  viscoelasticity, 
to  define  the  local  strain  energy.  This  approach  has  the  built-in  assump- 
tion that  there  is  a singular  relation  between  p and  the  boundary  displace- 
ments of  the  test  specimen. 

Figure  2 contains  photographs  of  a cut  In  a thin  piece  of  rubber 
dental  dam.  Figure  2a  shows  the  specimen  after  it  had  been  leaded  and  a 
notch  cut  into  It  (note:  the  crack  spontaneously  arrested  Itself  after 
propagating  a short  distance).  Figure  2b  is  a photograph  of  the  crack 
after  the  specimen  had  been  unloaded  and  reloaded  to  approximately  the 
same  load  level  as  seen  the  first  time.  Figure  3.  The  very  large  change 
In  the  radius  of  the  crack  tip,  p,  is  clearly  evident  from  the  two  photographs. 
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This  behavior  is  confirmed  for  solid  propellants  as  shown  in 
Figure  4.  This  test  was  performed  at  -65°F  and  Involves  a specimen 
that  was  repeatedly  loaded  to  cause  fatigue.  Figure  4a  shows  the 
virgin  specimen  at  the  point  of  failure  initiation.  Figure  4b  show.* 
the  same  conditions  for  the  fatigued  specimen. 

Andrews  (24)  was  the  first  to  observe  this  mixed  boundary  effect. 

Bills  (25)  made  the  additional  observation  that  the  rate  of  the 
second  loading  was  tne  primary  factor  determining  the  shape  of  the 
reloaded  crack.  Figure  5.  A low  rate  of  loading  is  requIrotTto  produce 
the  large  crack  tip  radius  noted  above,  see  Figure  5b.  But,  as  Bills 
observed,  the  unloaded,  pre-notched  rubber  specimen  has  a very  sharp 
crack  tip.  Figure  5a.  Thus,  a very  high  rate  of  loading  should  preserve 
the  sharp  crack;  which  it  does,  Figure  5c.  in  fact,  rapid  nulling  of 
the  sharp  crack  was  observed  to  reduce  rubber  specimen  strength  by  about 
a factor  of  10. 


Another  effect  of  strain  rate  upon  crack  tip  geometry  and  specimen 
strength  (In  a solid  propellant)  is  clearly  demonstrated  next. 

Effect  of  Strain  Rate  -a\  the  Crack  Tip 

The  mixed  boundary  effect  was  demonstrated  in  ANB-3066  propellant 
(CTPB).  The  double-notched  tensile  specisnan,  shown  in  Figure  6,  was  used; 
the  notches  being  razor  blade  cuts  iii  a single  pass  (to  prevent 
double-cutting  and  rounding  of  the  c. rack  tip}.  Five  different  notch 
depths  were  used:  0,  1/32,  T/16,  3/32  and  1/8  in. 


The  tests  were  conducted  at  a low  crosshead  rate,  so  that  the  very 
sharp  notch  effect  shown  In  Figure  4a  was  not  expected.  Instead,  reduced 
cross-section  should  experience  higher  strain  rates  with  a concomitant 
increase  in  the  apparent  strength  (based  on  the  net  specimen  cross-section). 
The  overall  specimen  had  a width  of  0.5  in.  Hence,  the  maximum  crack  depth 
of  1/8  In.  reduced  the  specimen  to  0.25  in.,  or  half  Its  original  cross- 
section.  All  of  the  failure  data  are  reported  in  terms  of  the  stress-at- 
rupture,  which  is  defined  In  terms  of  the  net  cross-sectional  area. 


The  local  strain  rate  at  the  crack  tip,  er»  may  be  crudely  estimated 
as  follows: 


, ,5b 


(2) 


Where  e,  is  the  overall  apparent  strain  rate  on  the  specimen 

a 

AQ  Is  the  original  cross-sectional  area  of  the  specimen 

A„  is  the  net  cross-section  at  the  notch 
n 
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a.  Initial  Tear  Contour,  Displacement  0.07  in. 


b.  Tear  Contour  after  Fatigue,  Displacement  0.30  in. 


Figure  4.  Tear  Contours  in  a Notched  Strip  Biaxial 
Specimen  of  Polyurethane  Propellant  at  -65°F 


a.  Unloaded  specimen 
with  cut  on  one 
side. 


b.  Shape  of  cut 
after  a low 
rate  of  loading. 


c.  Shape  of  tear 
tip  upon  rapid 
loading. 


figure  5.  Effect  of  Rate  of  Loading  Upon  Crack  Tip  Geometry 
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Thus,  we  estimate  that  the  notched  specimens  should  experience 
changes  In  strain  rate  up  to  a factor  of  about  2,  for  the  1/8  In.  notch 
depth.  At  77°  and  110°F,  the  strain  rate  effects  are  quite  small,  so  we 
expected  little  sensitivity  to  crack  depth.  Actually,  the  110°F  data 
gave  the  expected  results  (Figure  7),  while  the  77°F  tests  showed  more 
sensitivity  to  crack  depth  at  1/8  In.  than  was  predicted,  Figure  8. 

The  0°  and  -40°F  test  results  were  quite  sensitive  to  the  crack 
depth,  as  expected;  Figures  9 and  10.  Similarly,  at  -110°F  little 
sensitivity  to  strain  rate  was  expected  since  this  temperature  is  near 
the  glassy  temperature  of  the  propellant;  none  was  observed.  Figure  11. 

A large  statistical  variation  was  obtained  in  the  ~110°F  testing, 
however.  This  is  to  be  expected  since  glassy  materials  are  more  notch- 
sensitive;  and  our  methods  of  preparation  could  not  be  precise.  This 
becomes  abundantly  clear  when  one  considers  that  the  razor  blades  must 
drive  some  of  the  oxidizer  particles  ahead  of  it  and  into  the  cuts  making 
them  more  ragged,  and  effectively  deeper,  than  intended. 

Variations  In  the  test  data  may  also  arise  from  material  property 
gradients  in  the  carton  of  propellant,  from  which  the  specimens  were  cut. 

Figure  12  was  prepared  to  summarize  the  previous  data.  It  should 
be  clear  from  these  results  that  the  crack  merely  acts  like  a strain  rate 
magnifier,  as  previously  stated. 
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Crack  Depth,  in. 

Figure  8.  Strength  vs  Crack  Depth  in  ANB-3066  Propellant  at  77°F 

{Double-Notched  Specimens) 
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Figure  9.  Strength  vs  Crack  Depth  in  ANB-3066  Propellant  at  0°F 

(Double-Notched  Specimens) 
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Crack  Depth,  In. 

Figure  11.  Strength  v>  Crack  Depth  in  ANB-3G66  Propellant  at  -110°F 

(Double-Notched  Specimens) 
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Figure  12.  Strength  ys  Crack  Depth  in  ANB-3066  Propellant  at  Various 

Test  Temperatures 
(Double-Notched  Specimens) 
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SUPPORTING  GRIFFITH 

The  negative  Impressions  given  In  the  previous  section  must  be 
dispelled  somewhat,  since  there  are  conditions  where  a Griffith-type 
failure  criterion  seem  to  apply  to  propellant  failure  data. 

A completely  empirical  failure  criterion,  that  reduces  to  the 
Griffith  relation  when  the  propellant  Is  glassy,  is  described  next. 

Then,  using  that  relation  the  effects  of  crack  depth  in  propellant 
specimens  are  considered. 

Strength-Modulus  Relation  for  Creep-to- Failure  Data 

In  the  course  of  a number  of  Navy-Sponsored  programs  (7  to  11), 

Bills  demonstrated  that  propellants  seem  to  follow  a maximum 
principal  stress  failure  criterion.  According  to  this  approach,  the 
real  criterion  Is  the  tlme-to-fallure  of  a specimen  when  held  under  a 
constantly  applied  "true11  stress  (creep).  For  testing  conditions  other 
than  creep,  the  data  were  converted  to  the  "equivalent"  creep  condition 
by  means  of  a linear  cumulative  damage  relation  (26). 

A typical  failure  curve  Is  shown  in  Figure  13.  The  curve  was 
obtained  from  uniaxial  tests  perform  under  superimposed  hydrostatic 
pressures,  which  required  shifting  the  failure  data  along  the  time 
axis  by  a time-pressure  shift  factor,  ap.  As  shown  in  Reference  9, 
data  from  strip-biaxial  and  poker  chip  tests  fall  on  this  failure  curve. 

The  poker  chip  specimens  were  tested  at  Rocketdyne  and  contained  a 
center-mounted  stress  gage  that  easily  measured  the  local  stresses 
and  the  tlme-to-failure  initiation. 

The  curve  given  in  Figure  13  may  be  described  mathematically 
as  follows: 

o(?f)  - = fUf)  (4) 

Where  a{^)  is  the  applied  true  stress  expressed  as  a function  of 
is  the  reduced  time  to  failure  under  the  applied  stress 

o is  the  stress  below  which  failure  will  not  occur 

00 

fUf)  is  a descriptive  function  of  ^ that  need  not  be  defined  at  this  time 

The  reduced  time,  is  defined  by  the  relation 

Kf  = ^f/aTap  (5) 
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Figure  13.  Effect  of  Superimposed  Pressure  on  Constant  Stress  to  Failure  Data  for  a CTPB  Propellant 
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Where  ty  is  the  time  to  failure  under  a constant  stress 
ay  Is  the  time-temperature  shift  factor 

This  relation  may  be  normalized  with  respect  to  the  glassy  properties 
upon  division  by  og,  to  give: 

°(£f)  - 0OD 

— = gUf)  (6) 

°g  T 

Where  Og  Is  the  strength  In  the  glassy  region 
and 

g(«f)  = f(€f)/og  (7) 

The  relaxation  modulus,  E(c),  may  be  represented  by  the  following 
well  known  relation. 


E U)*Ee  +(Eg-Ee)_J  h(T)e*5/TdlnT  (8) 

Where  Eg  Is  the  equilibrium  relaxation  modulus 
Eg  is  the  glassy  modulus 

h(T ) is  a distribution  function  of  relaxation  times 
t is  a relaxation  time 


£ is  the  reduced  time  of  the  test 

After  noting  that  Eg»Ee,  Equation  (8)  can  be  rearranged  to  give  the 
following  normalized  relaxation  modulus 


EU)  - E£ 


w 

= J h(T)e"?/Tdinx 


(9) 


The  normalized  curves 
are  plotted  on  the  same  time 
to  be  made  are,  that  the  two 
region,  and  they  are  greatly 


for  relaxation  moduli  and  the  failure  stress 
scale  in  Figure  14.  The  most  notable  points 
curves  have  different  slopes  in  the  transition 
separated  along  the  time  coordinate. 
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thp  0f  ti]A  relaxat1on  modulus  data  and  shifting 

?!  ^^a!onVh  time-axis  permitted  the  fit  shown  In  Figure  15.  Here, 

of  t™  2at™dU,US  data  are  Seen  t0  fan  witMn  the  data  scatter 

equality^6  superpom1on  of  the  t"0  curves  in  Figure  15  permits  the  following 


aUf)  - ffa 


1/2 


(10) 


Where  E(cf/q)  Is  the  relaxation  modulus  taken  .jt  the  time  cf/q 

q I?"  efR?1r1cal  constant  that  shifts  the  entire  relaxation  modulus 

LUiVG 

The  value  of  q for  these  data  was  found  to  be  about  3 x 106.  This  is 
typical  of  past  data  where  values  of  10°  were  obtained. 

for  gl assySmaterfa?ngives*>e  Sll,ip,,f1cd  aftsr  "ot1"S  th=t  Griffith's  criteria 

1/2  1/2 


_ kr 


-1/2 


(11) 


Where  k Is  an  empirical  constant 

r is  the  free  surface  energy  per  unit  area  of  new  surface 
C is  the  length  of  an  internal  crack  In  the  specimen 

Inserting  Equation  (11)  into  (10)  and  rearranging  gives 

1/2  1/2 


- a-  ■ 5t? 


C«W  - M 


02) 


the  becoraes  a viscoelastic  fracture  criterion  that  includes 

* iriffith  re^a't^on  as  a special  case.  Furthermore,  the  Griffith  constants 
in  the  glassy  region  should  apply  over  the  entire  viscoelastic  time-scale. 

The  most  unique  feature  of  Equation  (12)  is  the  empirical  time-sHft 
A?"8  used  Anally  by  Halpin  (27)  and  given  a p^sic™  int'J- 
Sf  the  factor  m'  6 d0Ubt  tl,at  interPretat1on-  bu«  acknowledge  the  usefulness 
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Figure  15.  Superposition  of  Normalized  Moduli  and  Failure  Stress  Data 
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Recently,  Schapery  (28  and  29)  developed  a theoretical  failure 
theory  th^t  is  close  to  that  of  Equation  (12).  The  major  discrepancy 
is  in  the  exponent,  1/2,  on  the  modulus.  In  Schapery's  analysis  this 
exponent  has  a value  of  l/2(n+l);  where  n is  a constant.  Our  experi- 
mental data  support  the  possibility  that  Schapery's  relation  could  be 
correct. 

On  the  other  hand,  Blatz,  using  the  concepts  of  viscoelasticity, 
and  the  power  balance  concepts  of  Griffith  developed  a relation  that 
yields  a square  root  relation  (exponent  equal  to  1/2)  on  the  modulus. 

The  work  is  summarized  in  Appendix  A of  this  report. 

Tests  of  Crack  Depth  Dependence 

Two  tests  of  our  viscoelastic  failure  criterion  were  conducted. 

The  first  was  a rough  preliminary  effort  using  constant  load-to-failure 
tests  on  specimens  that  were  notched  on  one  side.  The  second  set  of 
tests  was  designed  to  eliminate  some  of  the  testing  errors  (specimen 
distortions  and  Inaccurate  time  measurements  obtained  In  the  first 
experiments.  As  shown  below,  both  experiments  qualitatively  support 
our  empirical  failure  criterion 

From  Equatif.i  (12)  It  Is  seen  that  a plot  of  [o(sf)  - oco]C^2  (13) 

versus  log  should  provide  a direct  means  for  evaluating  crack  depth 
dependence. 

Constant  Load-tp-Failure  Testing  of  Notched  Specimens 


Creep- to-fai lure  tests  were  made  at  77°F  on  specimens  with  razor 
cuts  of  0,  1/16,  1/8,  1/4  in.,  and  under  initial  stresses  of  45,  55,  65, 

90  and  100  psi,  based  on  the  net  cross-sectional  area  (after  the  cuts 
were  made). 

The  test  results  are  shown  in  Figure  16.  For  the  notched  specimens, 
the  ^ata  are  plotted  as  log  [oUf)  - aoJC’'2  versus  log  t^.  Although  there 
is  a considerable  scatter  in  the  data,  they  show  a general  consistency  with 
all  of  the  results  falling  along  the  same  curve. 
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Results  of  Constant  Load-to-Fallure  Tests  of 
Specimens  Notched  on  One  Side 
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The  notched  specimen  data  show  a strong  '.-pward  curvature  which 
is  mostly  associated  with  the  larger  stresses  and  the  deepest  cracks 
(1/4").  Therefore,  we  interpreted  the  curvature  to  be  an  artifact 
of  specimen  distortion  and  to  the  inaccuracies  of  measuring  failure 
times  below  10  sec. 

Results  of  the  unnotched  specimen  tests  are  plotted  as  log 
[o^f)  - a*,]  versus  log  tf.  They  follow  a straight  line,  as  expected. 

When  these  data  are  shifted  vertically  downward  until  they  superpose 
upon  those  of  the  notched  specimens,  we  can  obtain  an  estimate  for 
the  basic  flaw  size  in  the  unnotched  specimen.  The  dashed  line  in 
Figure  15  represents  the  vertically  shifted  curve.  The  crack  size 
represented  by  this  shift  is  about  0.032  in.  (800  u). 

Constant  Rate  Testing  of  Double-Notched  Specimens 

To  minimize  specimen  distortion  in  this  testing,  we  used  the 
double-notched  specimen  described  previously.  Figure  6.  To  minimize 
the  time  measurement  errors,  we  conducted  all  of  the  tests  under  a 
constant  displacement  rate  of  2 in. /min.  Also,  the  tests  were  con- 
ducted in  duplicate  at  -40°,  0°,  77°  and  110°F. 

The  notches  were  cut  into  both  sides  of  the  specimens  using 
single  passes  of  a razor  blade,  with  notch  depths  of  1/32,  1/16, 

3/32,  and  1/8  in.  In  this  planning  we  made  a mistake.  The  smallest 
notch  was  found  to  be  of  the  same  depth  as  the  naturally  occurring  flaw. 

At  the  same  time,  when  we  cut  into  the  specimen  we  force  oxidizer  particles 
to  "snow-plow"  ahead  of  the  razor  blade  effectively  cutting  deeper  than 
planned.  We  estimate  that  the  effective  crack  was  about  1/16"  deeper 
than  intended. 

The  experimental  data  were  reduced  to  give  "equivalent"  constant 
stress-to-failure  values  through  the  use  of  linear  cumulative  damage 
analyses  (3  and  6).  This  is  a correction  on  the  time-to-failure  that 
varied  from  0.15  to  0.58  for  these  results.  The  aT  values  were  taken 
from  past  data  on  AN8-3066  propellant. 

Figure  17a  contains  a plot  of  the  reduced  data  given  as  log 

[o(Sf)  _ ojc1/2  versus  log  tf/aT.  Also  included  In  Figure  17a  are  the 
failure  results  for  the  unnotchea  specimens.  The  slope  of  the  best 
straight  line  through  these  data  was  used  to  project  a line  through 
the  fracture  data. 

Comparing  the  actual  data  with  the  drawn  line,  in  Figure  17a, 
shows  marked  discrepancies  for  those  specimens  with  shallow  cracks. 

This,  we  feel,  is  due  to  the  specimen  preparation  (as  noted  above) 
and  to  the  relatively  equal  sizes  of  the  shallowest  crack  and  the 
inherent  flaws  (both  about  1/32  in.). 
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Unnotched  Specimens 
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Figure  17b  was  prepared  on  the  assumption  that  our  specimen 
preparation  actually  produced  cracks  that  were  1/16  in.  greater  than 
planned.  The  results  are  good,  seemingly,  justifying  our  preliminary 
assumption  of  the  flaw  size. 

In  spite  of  the  experimental  difficulties  the  initial  results 
certainly  give  support  to  the  frac*'_,-e  relation  (12),  and  justify 
further  investigations  of  the  -oar*  Admittedly,  these  tests 
should  be  pe/fortned  again:  , *er  circumstances  where  the  specimens 
and  the  cracks,  can  be  made  larger  to  eliminate  the  sensitivity  to 
specimen  preparations. 


Conclusions 

From  these  results  we  concluded: 

(1)  That  the  crack  depth  dependency  is  about  as  Griffith  had 
predicted,  even  under  viscoelastic  conditions. 

(2)  That  a constant  value  of  r may  characterize  propellant 
failures  over  the  entire  range  of  temperature  and  time  conditions 
where  propellants  are  normally  used.  This  is  not  expected  to  hold 
for  those  conditions  where  the  propellant  changes  its  physical  state 
(i.e.  crystallinity)  or  chemical  make-up  within  the  range  of  test 
temperatures. 

(3)  Reduction  of  the  failure  data  to  an  "equivalent"  creep- 
to-failure  test  may  take  into  account  crack  initiation  under  many 
complex  loading  histories.  Admittedly,  this  is  mainly  an  extension 
of  our  own  vast  experience  with  linear  cumulative  damage  testing. 

But,  once  we  have  clarified  the  role  of  crack  depth  dependence  the 
past  experience  becomes  totally  applicable. 

(4)  Future  work  should  be  directed  at  the  evaluations  of  r 
and  its  use  in  predicting  failures  under  various  loading  conditions. 

(5)  The  effects  of  "dull"  or  rounded  crack  tips  needs  to  be 
explored  and  corrections  included  in  the  basic  relation. 

(6)  The  strain-level  dependence  of  the  modulus  mav  be  a 
problem  when  evaluating  the  constancy  of  r in  Equation  (12).  This 
point  is  discussed  briefly  in  a later  section. 
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FACTORS  TO  BE  TAKEN  INTO  ACCOUNT  IN  DEVELOPING 
A FRACTURE  THEORY 

Developing  an  accurate  failure  criterion  for  solid  propellants 
will  be  a fantastically  complex  undertaking  since  a number  of  intricate 
processes  are  simultaneously  operating.  Our  recent  studies  have  clari- 
fied the  roles  of  a lumber  of  these  processes,  as  described  below. 

Basic  Failure  Ypchanisms 

Here  we  shall  consider  the  initiation  and  propagation  of  fractures 
in  solid  propellants  and  their  weak  binders. 

It  is  concluded  that  solid  propellants  will  initiate  fractures  in 
hydrostatic  tension,  but  propagation  requires  a shear  (distortion)  com- 
ponent. Furthermore,  propagation  of  the  cracks  is  such  as  to  make  them 
turn  and  stop,  necessitating  new  cracks  to  be  initiated.  The  consequence 
of  this  behavior  is  that  we  have  a "stick-slip"  situation  with  many 
failure  initiation  sites. 

These  features  are  described  next. 

Failure  in  Propellant  Binders 

The  failure  processes  In  a propellant  are  best  considered  in  terms 
of  the  binder  itself.  Consider,  first,  a simple  poker  chip  specimen  with 
a soft  transparent  binder  bonded  to  a glass  plate.  With  this  specimen 
we  can  view  failures  as  they  occur. 

Figure  18  is  a photograph  of  this  specimen  after  failure  initiation 
in  the  binder.  The  failures  are  seen  to  be  like  gas  bubbles  in  a liquid. 

Tbe  bubbles  do  not  seem  to  propagate  significantly  after  their  initial 
appearance  (although  they  may  elongate  somewhat  as  the  specimen  is  stretched). 
Instead,  there  seems  to  be  a tendency  to  initiate  new  failure  sites. 

Propagation  of  these  failures  clearly  requires  an  additional  driving 
factor.  We  found  this  in  a test  where  a poker  chip  specimen  experienced 
combined  tension  and  shear.  As  shown  by  the  photograph  in  Figure  19  the 
bubbles  have  grown  in  the  direction  of  the  shear  component.  The  same 
failure  result  was  obtained  in  a clear  liner  that  was  bonded  to  a glass- 
encased  5 in.  dla.  grain,  Figure  20. 

When  the  bubbles  form  under  conditions  of  combined  tension  and  shear, 
there  is  a peeling  action  in  the  direction  of  the  shear.  In  the  motor  this 
is  seen  as  long  tubular  separations  that  wander  In  fern-like  patterns. 

Figures  21  and  22.  The  fern-like  pattern  will  continue  to  grow,  since  the 
end  stress  condition  (combined  radial  tension  and  axial  shear)  moves  along 
with  the  leading  edge  of  the  growing  crack. 


Preceding  page  blank 
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This  mechanism  of  crack  propagation  Is  directly  relevant  to 
propellant  failure  processes.  The  close  packing  of  oxidizer  particles. 
Figure  23,  produces  local  poker  chip-type  conditions,  which  Leon  and 
McClintock  (30)  have  show”  to  produce  stress  levels  up  to  twice  those 
imposed  on  the  overall  specimen.  Failure  is  induced  in  the  propellant 
at  low  elongations  through  bubble  formation  in  the  binder  or  through 
failure  of  the  oxidizer-binder  bond  (usually  called  dewetting  in  the 
industry).  This  initial  failure  is  not  catastrophic  to  the  propellant, 
so  it  will  continue  to  elongate  expanding  the  vacuoles  previously  formed. 

If  we  concentrate  on  the  dewetting  behavior  at  a crack  tip  the 
vacuoles  become  very  highly  elongated,  as  shown  in  Figure  24.  At  the 
same  time,  the  contiguous  vacuoles  experience  strong  shearing  forces 
due  to  the  geometries  of  their  boundaries.  As  stated  previously,  this 
shearing  action  propagates  the  cracking  process  between  the  vacuoles. 

An  artifact  of  this  behavior  is  seen  through  the  microscope,  where 
these  ligaments  appear  to  be  flowing  (shearing  action)  like  a viscous 
liquid. 


Localized  vacuole  formation  has  many  implications  in  the  propellant 
failure  process,  including  a marked  effect  upon  the  radius  of  curvature 
of  the  crack  tip.  This  is  discussed  further  in  a later  section. 

Appearance  of  Real  Cracks 

Two  distinct  failure  processes  may  be  defined  for  solid  propellants. 
They  are  defined  by  their  propagation  rates,  with  the  low  rate  processes 
giving  jagged  cracks,  while  those  at  very  high  rates  yield  brittle-type 
failures. 

All  observations  of  inner-bore  grain  cracking  under  low  deformation 
rates,  have  shown  multiple  tear  initiations.  This  is  easily  seen  in  the 
early  stages  of  these  failures  when  the  inner-bore  frequently  develops 
minute  cracks.  As  shown  in  Figures  25  and  26  these  tiny  cracks  propagate 
a short  distance,  then  turn  from  their  original  trajectories  and  stop. 

The  final  joining  of  these  short  tears  yields  the  jagged  cracks  shown. 

(1)  Low  Propagation  Rates 


All  observations  of  inner-bore  grain  cracking  have  shown  multiple 
tear  initiations.  This  is  easily  seen  in  early  stages  when  the  inner-bore 
frequently  develops  minute  cracks.  As  shown  in  Figures  24  and  25  these 
tiny  cracks  propagate  a short  distance,  then  turn  from  their  original 
trajectories  and  stop.  The  final  joining  of  these  shon;  tears  yields 
the  jagged  cracks  shown. 

Non-linear  propagation  of  grain  cracks  also  occur  in  the  radial 
direction,  Figure27.  Here,  a 6 in.  diameter,  cyllndrically  perforated, 
grain  was  failed  by  internal  pressurization.  Figure  27  is  a photograph 
of  the  central  cross-section  of  the  grain. 
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F*gure  25. 
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Figure  27.  Non-Linear  Propagation  of  Radial  Cracks 
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Other  examples  of  multiple  fracture  Initiations  and  jagged  cracks 
are  given  In  References 

A partial  explanation  for  this  behavior  may  be  obtained  from  photo- 
elastic analyses.  Figure  28  shows  the  photoelastic  stress  patterns  at  the 
crack  front  In  a rubber  binder.  Figure  28a  was  obtained  with  the  stress 
at  a sufficiently  low  level  that  the  crack  remained  stationary,  while 
Figure  28b  was  made  at  a level  where  the  crack  was  slowly  progressing. 

The  artificial  crack  was  first  Introduced  by  a sharp  razor  cut  at  the 
edge  of  the  thin  specimen.  It  is  seen  that  the  sharp  front  acquires  a 
significant  radius  of  curvature  before  the  front  progresses. 

On  close  examination  the  photoelastic  stress  patterns  show  three 
stress  concentration  points,  one  directly  in  front  of  the  crack  tip, 
and  two  near  the  tip,  but  at  an  angle  to  the  direction  of  crack  propaga- 
tion. 


Using  chese  findings,  and  a modification  of  a concept  by  Andrews 
(31  and  32),  we  have  concluded  that  the  stress  at  the  axis  of  the  crack 
tip  is  attenuated  by  the  dewetting  (or  vacuole)  process,  leaving  the  two 
lateral  stress  concentrations  to  fight  for  control  of  the  growing  crack. 
Thus,  crack  propagation  will  not  follow  the  crack  axis  but  will  follow 
one  of  the  maximum  stress  trajectories,  see  Figure  29.  This  is  what 
appears  to  happen  in  the  small-scale  propagation  of  a crack  In  propellants. 

In  the  grain  the  rotation  of  the  crack  turns  it  into  the  hoop 
direction  where  the  normal  stress  conditions  are  insufficient  to  continue 
the  propagation  process.  Thus,  the  crack  terminates  and  a new  crack  must 
be  initiated  to  continue  the  failure  process. 

On  the  basis  of  the  above,  the  cracking  of  a propellant  grain  is 
seen  to  go  through  four  stages:  initiation,  propagation,  and  termination 
of  small  crack  segments  followed  by  a final  tearing  process  across  the 
various  small  cracks  to  form  one  large  crack.  Acini ttedly,  stress  levels, 
test  temperatures  and  material  properties  will  modify  the  various  steps 
leading  to  a wide  variety  of  behaviors  from  smooth  to  very  jagged  crack 
surfaces. 

(2)  High  Propagation  Rates 

Under  special  condi t’*ons  of  temperature,  deformation  rate,  or, 
sometimes,  depth  of  crack,  toe  propellant  will  exhibit  a high  rate  of 
crack  propagation.  Figure  30  is  a photograph  of  an  inner-bore  crack 
that  was  generated  at  a moderate  rate.  The  crack  is  seen  to  be  slightly 
rough,  but  it  is  not  jagged  like  those  obtained  at  lower  rates. 

Under  brittle  failure  conditions  where  the  crack  rates  are  very 
high,  the  failures  are  smooth  and  they  propagate  through  the  AP  crystals, 
as  well.  We  have  no  available  photographs  of  this  phenomenon.  But,  a 
bar  of  propellant  soaked  in  liquid  nitrogen  will  exhibit  brittle  fracture, 
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Axial  Direction  of  the 


of  Actual  Crack 


Figure  29.  Deviation  of  Growing  Crack 
from  Axis  of  Symmetry 
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Figure  30.  Long  Crack  Formed  at  an  Intermediate 
Propagation  Rate 
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and  the  test  Is  readily  conducted  by  hand  (gloved,  of  course).  Actually, 
most  propellants  propagate  fractures  so  readily  at  this  temperature  that 
a JANNAF  tensile  bar  can  be  reduced  to  powder  upon  crushing  in  the  fingers 
of  one  hand. 

(3)  Conclusions 

The  various  failure  processes  observed  above  lead  to  the  following 
conclusions: 

(a)  Under  low  propagation  rates  there  will  be  a considerable  energy 
dissipation  along  the  jagged  crack  front,  and  in  the  extraneous  failure 
initiations  away  from  the  main  crack. 

(b)  The  non-linear  propagation  of  cracks  at  low  rates  should  lead 
to  errors  in  predicting  crack  trajectories. 

(c)  At  very  high  propagation  rates  the  crack  trajectories  are  linear 
and  the  failed  surfaces  are  smooth. 

(d)  In  developing  a practical  theory,  the  transition  from  low  to 
high  propagation  rates  must  be  accomplished  by  analytical  descriptors  of 
the  associated  mechanical  phenomena. 


Effect  of  Dewetting  Upon  the  Criterion  of  Failure 

For  specimens  not  experiencing  compressive  strains  the  dewetting 
phenomenon  (or  vacuole  formation)  provides  a simple  control  over  the  mode 
of  failure.  In  effect,  dewetting  reduces  the  local  ligament  geometry  to 
that  for  an  equivalent  uniaxial  test,  regardless  of  the  loading  history. 
There  is  support  for  the  same  behavior  under  lateral  compression,  but  it 
will  not  be  made  at  this  time. 

Uniaxial,  Biaxial  and  Trlaxial  Tens ion 


Figure 31  is  a photomicrograph  of  dewetting  cavities  (dark  bands) 
forming  around  some  glass  bead  as  the  specimen  is  uniaxially  strained. 

Figure  32  is  another  photomi crograph  of  a specimen  that  has  been  pulled 
equally  in  biaxial  tension.  (This  was  an  "Iron  Cross"  specimen  that  shows 
extensive  ligamentati on  in  a portion  of  a leg  that  was  uniaxially  strained.) 
As  shown  in  Figure  33,  the  dewetting  forms  a honey-comb  pattern.  Thus, 
at  the  local  level,  the  dewetting  is  the  same  as  that  for  uniaxial  tension. 

This,  reduction  of  problems  to  simple  uniaxial  dewetting  conditions 
helps  to  explain  why  the  uniaxial,  biaxial  and  poker  chip  data  can  be 
represented  by  the  same  Maximum  Principal  Stress  failure  criterion  (see 
above). 
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DIRECTION  OF  STRAIN 


Figure  33.  Schematic  Representation  of  Dewetting  Bands  in  a Biaxially  Strained  Specimen 


Shear 


NWC  TP  5684 


We  ran  a special  test  on  a shear  specimen  using  notches  cut  deeply 
into  both  ends.  Figure  34.  Extensive  dewetting  occurred,  especially  at 
the  right  h.nd  notch.  The  dewetting  bands  fell  at  an  angle  to  an  imaginary 
line  connecting  the  notches.  Thus,  even  the  shear  specimen  yields  extensive 
dewetting;  the  behavior  being  like  a rotation  of  the  uniaxial  tensile  test. 

Conclusion 


For  the  usual  conditions  of  solid  propellant  applications,  there 
is  extensive  dewetting  with  the  result  that  all  testing  modes  act  as  if 
they  were  uniaxials.  However,  at  low  temperatures  dewetting  may  be  sup- 
pressed and  the  simplifying  assumptions  due  to  the  dewetting  bands  will 
no  longer  hold. 


Loading  Path  Dependence 

Failures  in  solid  propellants  depend  greatly  upon  the  order  of 
application  of  the  loads.  Examples  of  this  behavior  are  plentiful. 
Those  presented  here  were  selected  to  represent  path-dependence  under 
several  different  types  of  loading. 

Temperature-Strain  Paths 

In  this  example  we  have  considered  three  paths:  Path  I is  for  a 

propellant  specimen  that  is  strained,  then  cooled;  Path  II  is  for  the 
specimen  that  is  cooled  and  then  strained;  and  Path  III  Involves  simul- 
taneous cooling  and  straining.  Figure  35  summarizes  these  paths 
schematically. 


The  path  dependence  of  a polyurethane  propellant  was  observed  in 
a series  of  tests  where  the  specimen  was  strained  or  cooled  along  the 
different  paths  of  Figure  35.  The  results  are: 

Path 

Description 

Observed  Failure 
Strains 

I 

Cool  specimen  from  +40°  to 
-75°F,  then  strain  to  break 
at  0.74  min"' 

Average  e.  = 2.0% 
at  -75°F  D 

II 

Stretch  specimens  to  various 
strain  levels  at  40°F  hold 
and  cool  to  -75°F 

Strains  up  to  25% 
held  9 days  when 
test  terminated 

III 

Simultaneous  cooling  and 
straining  from  +40° F to 
-75°F  at  10°F/hr  and  1.08% 
strain/hr. 

Average  ck  = 7.7% 
at  -30°F 

III 

(Modified) 

Specimens  were  given  a 10% 
pre-strain  before  following 
Path  III  above. 

Average  eb  = 21.9% 
at  -75°F 

-55- 


a.  Unstrained 


b.  20%  Shear  Strain 

Figure  34.  Dewetting  in  a Notched  Specimen  Undergoing  Shear  Deformation 
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From  these  data  it  can  be  seen  that  Path  I gives  the  lowest  value 
for  cfc,  while  Path  II  gives  the  largest  values.  However,  Path  III  is  more 
typical  of  the  experience  of  a motor,  and  is  Intermediate  between  Paths  I 
ard  II. 

Bauschinger  Effect 

This  effect  modifies  the  response  behaviors  of  all  materials  that 
exhibit  creep  or  hysteresis,  including  hard  steels.  The  effect  described 
here  arises  from  a pre-loading  in  compression  during  which  the  material 
experiences  "temporary  set."  This  "set"  causes  an  effective  shift  of  the 
coordinate  stress -strain  axes,  which,  of  course,  will  modify  the  apparent 
failure  properties  wnen  the  sample  is  tested  in  tension. 

Tests  were  conducted  on  the  standard  specimen  used  to  characterize 
the  strain  energy  release  rate  in  solid  propellants.  This  is  a strip- 
biaxial  specimen  with  a 1 in.  razor  cut  made  in  the  specimen  center. 

In  this  test,  the  specimen  was  put  into  5%  compression,  held  10  min. 
to  induce  a temporary  set,  then,  without  removing  it  from  the  fixture,  the 
specimen  was  pulled  at  a constant  rate  until  crack  initiation  was  detected. 

Figure  36  contains  a plot  of  the  stress-displacement  curve  for  a 
CTPB  propellant  tested  in  this  manner.  The  results  of  a standard  test 
are  included  for  comparison  in  Figure  36.  As  summarized  in  Table  I there 
Is  a reduction  of  15%  in  the  failure  stress  (tear  initiation)  while  the 
axial  deformation  is  reduced  by  24%. 

Pre-Strain  Before  Pressurization-to-Failure 

This,  also,  is  a Bauschinger  effect,  except  that  the  shift  of  the 
coordinate  axes  is  into  the  tensile  stress -strain  region.  Thus,  it  may 
show  a benefit  to  grain  performance. 

Motor  firings  at  low  temperatures  involve  rapid  straining  beginning 
from  the  existing  thermal  strain.  The  laboratory  test  for  this  is  to  pre- 
strain the  propellant  then  test  it  to  failure  under  superimposed  hydrostatic 
pressures. 

^9u|,,e  37  contains  the  failure  envelopes  for  a CTPB  propellant  tested 
at  -75  F under  an  800  psl  superimposed  hydrostatic  pressure.  The  envelope 
generated  from  specimens  with  a 10%  pre-strain  has  a distinctly  improved 
failure  strain  capability  over  that  of  the  specimens  without  pre-strain. 

Figure  37  was  used  to  plan  a set  of  experiments  where  four  motors 
were  cooled  to  -75°F  and  pressure  tested  (33).  Failures  were  not  obtained 
until  the  envelope  on  the  right  was  exceeded. 
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rioii-f‘37.  Effect  of  Pre-Stnin  on  the  Failure  Envelope 
+'o>'  a Crr’!-  1 1 ant  Tested  Under  800  psi  Pressure  at  -75°F 
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Conclusions 


The  strong  path  dependency  of  solid  propellant  failures  is  primarily 
a material  response  problem.  These  changes  in  material  response  affect  both 
the  shape  and  the  stresses  at  the  tip  of  a crack,  thus  modifying  propellant 
failure  properties. 


If  we  are  to  have  a viable  failure  theory  it  must  contain  a realistic 
constitutive  theory  that  accounts  for  path-dependencies  like  those  described 
above.  Because  of  that  need  we  instituted  a theoretical  effort  to  develop 
the  required  relations.  The  preliminary  efforts  in  ibis  area  are  summarized 
in  a later  section  of  this  report. 

Independence  of  Filler  Content 

The  important  aspects  of  dewetting  and  vacuole  formation  need  further 
clarification.  This  is  essential  since  a number  of  failure  theories  (15  to 
17)  have  centered  around  the  phenomenon.  We  find  the  following  data  to  be 
somewhat  challenging  to  those  theories. 

A number  of  HTPB  propellant  compositions  were  made  using  different 
filler  fractions  (0,  50,  64  and  76.3  vol.  %)  and  various  binder  equivalents 
ratios  (crosslinking  levels).  One  of  the  material  characterization  methods 
included  standard  tensile  tests  conducted  over  a range  of  temperatures  and 
crosshead  speeds. 

The  tensile  strength  data  given  in  Figures  38  to  40  show  this  pro- 
perty to  be  entirely  independent  of  filler  content.  This  contradicts  our 
pre-conceptions  of  dewetting,  where  material  strength  would  seemingly  be 
reduced  as  the  filler  content  increases. 

That  is,  one  would  expect  that  as  more  filler  is  used  there  should 
be  more  dewetting;  or,  that  the  particles  should  be  brought  closer  together 
inducing  local  stress  concentrations  and  earlier  failures.  Clearly,  the 
data  in  Figures  38  to  40  contradict  that  expectation. 

Strong  Time  Dependence 

There  are  a number  of  meaningful  molecular  parameters  that  affect  the 
way  propellants  respond  and  fail.  We  shall  consider  only  one,  an  old  "friend" 
binder  crosslinking,  which  we  have  found  to  be  much  less  important  than  pre- 
viously supposed. 

An  exception  is  taken  to  the  position  of  Landel  and  Fedors  (34  and  35) 
with  respect  to  the  failure  properties  of  lightly  crossl inked  elastomers, 
which  is  the  case  for  solid  propellant  binders.  They  indicate  that  the 
tensile  strength  of  an  elastomer  is  proportional  to  its  crosslink  density. 
However,  from  the  data  of  Bueche  and  Dudek  (18,  19  and  20)  this  should  be 
true  only  for  the  more  highly  crosslinked  elastomers. 

The  analyses  of  Bueche  and  Dudek  (18  to  20)  show  that  the  stress  at 
break  on  the  rubber  consists  of  two  parts;  a normal  kinetic  theory  stress  on 
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Master  Plot  of  the  Stress-at-Break  for  an  HTPB  Polymer 
- Equivalents  Ratio  60/40  - 


39.  Master  Plots  of  the  Stress-at-Break.  for  an  HTPB- Polymer 


Figure  40.  Master  Plots  of  the  Stress-at-Break  for  an  HTPB  Polymer 

- Equivalents  Ratio  80/20  - 
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the  network  and  a highly  rate- dependent  stress  held  by  fully  extended 
chains.  The  fraction  of  the  chains  in  the  fully  extended  state  at  the 
instant  of  break  is  much  greater  for  lightly  crosslinked  rubber  than  for 
highly  crosslinked  rubbers  under  nonequilibrium  test  conditions.  They 
show  that  the  kinetic  theory  stress  (related  to  polymeric  crosslinking) 
contributes  the  major  portion  of  the  stress-at-break  in  highly  crosslinked 
amorphous  rubbers.  However,  the  stress-at-break  of  very  lightly  cross- 
linked rubbers  is  largely  the  result  of  stresses  in  the  most  highly 
elongated  chains. 

This  theory  is  illustrated  in  Figure  41,  which  is  a plot  of  the 
stress-at-break  vs  the  log  of  the  crosslink  density  for  an  SBR  amorphous 
rubber.  The  experimental  data  show  the  strength  to  increase  to  a maximum 
then  decrease  as  the  crosslink  density  is  increased.  The  theoretical 
prediction  of  Bueche  and  Dudek  is  given  as  the  solid  line  in  Figure  41 
while  the  load  supported  by  the  crosslinked  network  (elastic  contribution) 
is  shown  by  the  dashed  line. 

The  maximum  in  the  curve  of  Figure  41  was  also  studied  by  Bueche  and 
Dudek  (20)  and  an  explanation  given.  It  is  due  to  the  fact  that  a greater 
fraction  of  chains,  which  have  been  extended  to  their  limit  of  extensi- 
bility, is  present  at  the  moment  of  fracture  in  lightly  crosslinked  vulcani- 
zates.  A rate  dependent  nonaffine  deformation  of  the  network  junction 
points  comes  into  play  at  these  chain  elongations  approaching  full  exten- 
sion and  provides  a molecular  mechanism  to  explain  the  results. 

From  these  data  we  concluded  that  the  parameter  plays  a lesser  role 
in  propellant  properties  and  aging  effects  than  previously  supposed. 

Characterization  of  Material  Fracture 

Laboratory  characterizations  of  propellant  fracture  are  considered  to 
be  inaccurate.  This  stems,  mostly,  from  differences  between  the  shapes  of 
the  crack  tips  as  seen  in  the  motor  and  those  obtained  during  laboratory 
characterization. 

All  of  the  cracks  obtained  in  motor  testing  have  sharp  tips,  Figures 
25  to  27  and  Figure  30.  By  contrast,  the  crack  tip  in  a laboratory  specimen 
can  be  markedly  rounded,  Figure  42.  The  photograph  given  in  Figure  42  was 
was  selected  as  the  most  exaggerated  case  of  those  available  to  us. 

Figure  43  is  a photograph  of  fracture  initiation  in  a specimen  where 
the  razor  cut  was  slightly  canted. 

Figure  44  was  chosen  to  illustrate  that  the  stresses  at  the  crack  tip 
are  not  necessarily  the  highest  ones  in  the  test  specimen.  Here,  specimen 
failure  initiated  at  a point  ahead  of  the  artificial  crack  tip. 
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Figure  41.  Strength  of  Unextracted  SBR  at  an  Extension 
Rate  of  5 min.^.  (Copied  fron  Bueche  and  Dudek). 
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Figure  43.  Shape  of  Crack  Where  the  Initial  Razor  Cut  was  Slightly  Canted 

(Polyurethane  Propellant  at  -65°F) 
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Figure  44.  Failure  Initiation  Ahead  of  the  Crack  Tip  In  a Polyurethane 

Propellant  at  -65°F 
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In  summary,  we  feel  that  a material  characterization  method  that 
gives  sharp  cracks  must  he  used  in  the  future. 

Dependence  of  Low-Strain  Moduli  on  Dewetting 

Dewetting  in  propellants  appears  very  early  in  the  straining  process. 
This  is  readily  observed  in  careful  stress-strain  measurements  of  end-bonded 
test  specimens.  A typical  st.--,ss-strain  curve  is  illustrated  in  Figure  45. 

It  is  concave  downward,  the  modulus  continually  decreasing  with  increasing 
strain.  Schapery  (36)  noted  this  behavior  of  propellant  and  made  tests 
showing  that  stress  relaxation  moduli  taken  at  0.25%  strain  were  four  to 
six  times  those  taken  at  3%  strain. 

Francis  (21)  separately  noted  the  effect  of  strain  level  upon  stress 
relaxation  moduli,  while  Adicoff  and  Lepie  (22)  drew  the  conclusion  that 
the  dewetting  effect  causes  the  propellant  to  deviate  from  thermorheological 
simplicity.  In  particular,  they  noted  that  the  time-temperature  shift  factor, 
ay,  is  modified  by  the  strain  level. 

With  respect  to  Equation  (1)  this  strain  level  dependency  is  of  major 
importance.  If  r is  to  be  evaluated  as  a valid  constant,  it  will  be  necessary 
to  obtain  accurate  values  of  the  relaxation  moduli. 

It  is  recommended  that  the  various  available  techniques  be  evaluated 
for  accurate  moduli  determinations  at  strain  levels  in  the  vicinity  of  0.25% 
strain. 
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THEORY  - INTRODUCTION 

(A  New  Approach  to  Solid  Propellant  Mechanical  Behavior) 

Basic  Response  Processes 

Solid  propellant  material  is  a composite  of  rigid  incompressible 
filler  particles  (roughly  spherical)  entrained  in  an  Incompressible  rubber 
matrix.  When  such  a composite  is  deformed  mechanically,  the  rubber  stores 
elastic  strain  energy  and  the  binder  pulls  away  from  the  filler,  producing 
a modest  amount  of  dilatation  In  the  composite  structure.  In  addition  the 
rubber  may  dissipate  the  stored  energy  viscoelastically  and  the  filler 
particles  may  slip  relative  to  the  binder,  causing  additional  dissipation. 

In  the  theory  that  we  present  below,  these  dissipation  effects  are  neglected; 
the  intent  is  to  extend  the  theory  in  a later  development  program  to  include 
such  effects.  Thus,  we  are  talking  about  rigid  dewettable*  filler  particles 
in  an  incompressible  rubber  matrix. 

We  had  available  simple  tensile  data  (all  at  77°F)  at  various  strain 
rates  from  0.0125  to  1.25  min“1  and  at  various  superimposed  pressures  from 
0 to  1000  psig.  In  the  following  section,  we  shall  present  a constitutive 
theory  which  correlates  these  data  and  determines  the  five  parameters  of  the 
theory  in  a self-consistent  fashion.  The  theory  is  relatively  simple  and 
completely  physically  motivated.  In  a nt. Isbell,  it  postulates  that  the  pul  1 - 
away  of  the  binder  from  the  filler  is  controlled  by  the  Griffith  surface 
energy  criterion. 

Technical  Background 

The  science  of  solid  propellant  mechanic?-  behavior  has  had  a long  and 
fruitful  history.  In  1959,  workers  at  Aerojet-General  SRP  painted  circules 
on  simple  tensile  specimens  and  photographed  the  subsequent  deformation  of 
these  circles  onto  elipses  during  simple  tensile  stretch.  It  was  determined 
that  Poisson's  ratio  decreased  from  an  initial  value  of  1/2  to  a value  of 
about  1/4  at  40-50  percent  strain.  It  was  roughly  in  this  period  of  time 
that  workers  in  the  field  began  to  realize  that  the  fall  off  in  Poisson's 
ratio  and  the  concatenate  softening  of  the  modulus  (slope)  of  the  tensile 
curve  were  in  some  way  related  to  the  pull away  of  the  binder  from  the  filler, 
or  to  dewetting,  as  it  is  commonly  called. 

Roughly  ten  years  later  there  appeared  on  the  scene  another  Aerojet 
worker  (Farris)  who  proposed  a simple  constitutive  equation  to  explain  this 
behavior  (37,  33)  namely: 

0 = Ee(l  - 8~  ) (13) 

e 


Although  dewetting  is  used  for  simplicity  In  the  descriptions,  it  should 
be  apparent  that  tearing  of  *he  binder  adjacent  to  the  filler  will  lead 
to  the  same  descriptive  processes,  and  the  same  forms  for  the  relations 
generated  below. 
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This  was  the  first  example  of  a nonlinear  elastic  constitutive  equation 
based  on  ideas  other  than  those  whi  h stem  from  a strain  energy  function 
(it  Is  easy  to  show  that  Equation  (13)  is  not  analytic).  The  importance 
of  introducing  a non-analytic  function  of  this  sort  is  that  it  made  people 
realize  that  some  other  type  of  mathematics  (besides  finite  elasticity) 
was  needed  to  represent  propellant  behavior.  Fitzgerald  later  suggested 
(1968)  that  one  use  semi-norms  to  handle  the  representation.  Farris  (39) 
working  for  Fitzgerald.,  put  the  idea  on  a physical  basis  by  suggesting 
that  the  applied  strain  in  a tensile  test  be  constructed  in  a distributed 
fashion,  and  that  the  distribution  be  limited  by  the  ratio  of  a critical 
strain  to  a semi-norm  of  the  history  of  the  strain.  Farris'  rationali- 
zation of  the  use  of  the  semi-norm  of  the  history  of  the  strain  is  that 
it  is  a chronological  counter  of  damage. 

Taking  this  cue  from  work  such  as  Farris',  Lindsey  (40)  suggested 
that  an  alternative  approach  to  representing  the  damage  mechanism  existed 
in  the  form  of  classical  plasticity.  An  objection  exists  to  the  use  of 
plasticity  however  from  the  standpoint  that  a propellant  can  evidence 
marked  hysteresis  and  no  permanent  set. 

Farris'  idea  of  a chronological  counter  of  damage  in  the  form  of 
a semi-norm  has  proved  very  fruitful.  This  same  idea  however  has  led  to 
the  introduction  of  a large  number  of  semi -norm  terms  of  various  orders 
which  do  not  easily  line  up  with  physical  concepts.  It  is  the  desire  to 
provide  a completely  physical  theory  that  leads  the  present  Investigator 
to  propose  a damage  criterion  based  on  completely  physical  ideas  and 
physi cal ly-moti vated  mathematics . 

Furthermore,  it  is  recognized  by  many  workers  in  the  field  today 
that  it  Is  possible  to  replace  classical  plasticity,  with  Its  yield 
criterion,  by  a completely  differential  theory  which  involves  a semi- 
norm in  the  form  of  a chronological  counter.  For  example,  in  the  case 
of  simple  tension,  following  workers  like  Rivlin,  Pipkin,  and  Farris, 
one  can  write 

0 t 6 |e  jo  = Eg  (14) 

and  generate  hysteresis  with  no  discontinuity  in  stress  during  loading. 
Equation  (14)  is  just  one  example  of  many  forms  of  the  theory  developed 
by  Pipkin  and  Dag  (41)  using  chronological  variables  of  the  form: 

s =V?d~ 

Thus  it  appears  that  whether  one  Introduces  the  damage  mechanism  via 
plasticity  or  via  a Mull  In's  effect,  one  must  use  chronological  variables 
which  behave  like  semi -norms. 
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In  the  theory  that  we  are  about  to  propose,  the  role  of  the  semi- 
norm  Is  played  by  a physical  parameter,  namely  the  porosity  that  is  generated 
by  the  pull  away  of  the  binder  from  the  filler.  In  this  process  of  pullaway, 
new  surface  area  is  created  and  new  surface  energy  is  added  to  the  deformed 
propellant.  We  use  the  Griffith  criterion  to  control  the  rate  of  increase 
of  the  new  surface  energy. 

In  this  approach  to  fitting  data  based  on  our  theory,  we  shall  assume 
that  all  the  hysteresis  is  generated  by  the  irreversibility  of  the  porosity. 
We  shall,  therefore,  use  no  hereditary  integral,  but  will  use  a nonlinear 
measure  of  strain. 
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ANALYTICAL  DEVELOPMENT 
Generating  the  Physical  Model 

Roughly  speaking,  the  maximum  strain  evinced  by  a solid  propellant 
specimen  subjected  to  simple  tension  is  about  50  percent.  Within  this 
range  of  strain,  only  a mildly  nonlinear  constitutive  equation  Is  needed 
to  fit  the  behavior  of  a purely  elastic  material.  This  statement  implies 
the  presence  in  the  constitutive  equation  of  only  two  moduli,  the  shear 
modulus  and  the  bulk  modulus.  Because  the  binder  pulls  away  from  the 
filler  particle,  these  two  moduli  soften  with  increasing  strain.  This 
demands  that  the  theory  be  even  more  nonlinear.  One  way  of  handling  this, 
as  suggested  by  Farris  and  Schapery  is  to  make  the  shear  modulus  a 
decreasing  function  of  the  dilatation,  and  the  bulk  modulus  a decreasing 
function  of  the  octahedral  shear  strain.  Such  an  approach  is  subject  to 
two  limitations.  One  is  that  the  presently  proposed  forms  such  as  Equation  (13) 
are  not  consistent  with  linear  theory  in  the  limit  of  infinitesimal  strain, 
while  other  forms,  which  are  based  on  functions  that  are  essentially  empirical, 
lack  physical  motivation. 

In  an  attempt  to  motivate  the  decrease  of  the  shear  modulus  and  the 
bulk  modulus  by  a physical  argument,  we  propose  the  following  model. 

We  envisage  the  propellant  rubber-filler  matrix  as  a collection  of 
unit  spherical  cells  each  with  a dewetting  cavity.  As  the  propellant  Is 
strained,  the  cavity  grows  in  proportion  to  the  amount  of  new  surface  area 
that  Is  produced  when  the  binder  pulls  away  from  the  filler  (or  tears  in 
the  rubber  adjacent  to  the  filler).  In  effect,  our  model  treats  the  binder- 
filler  matrix  as  a porous  material  in  which  porosity  Increases  during  loading, 
remains  fixed  during  unloading  (this  makes  it  a physical  semi-norm),  and  in 
which  the  porosity  effects  the  bulk  modulus  and  shear  modulus  In  such  a way 
that  they  decrease  with  increasing  porosity.  In  turn,  the  rate  of  increase 
of  the  porosity  with  Increasing  strain  is  controlled  by  the  Griffith  criterion. 

The  physical  model  described  above  led  to  a straightforward  mathematical 
relation.  Testing  of  that  model  with  actual  propellant  response  data  required 
only  a minor  correction  for  the  case  of  superimposed  hydrostatic  pressures. 
These  analytical  relations  are  described  below. 

Choosing  the  Strain  Energy  Function 

In  dealing  with  nonlinear  materials,  there  are  a plethora  of  strain 
energy  functions  that  can  be  used  to  describe  elastic  behavior.  This  is 
particularly  true  In  dealing  with  materials  (such  as  propellants)  in  which 
the  tensile  strain  does  not  greatly  exceed  50%.  This  nonlinearity  is  mild, 
and  many  functions  describe  elastic  behavior  up  to  50%  with  equal  accuracy. 

For  the  initial  studies  we  chose  the  following  function: 


Preceding  page  blank 
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W = AI  + BA  (16) 

c 

Where  W is  the  elastic  strain  energy  function 
I is  the  first  strain  invariant 

G 

A is  the  dilatation 
A and  B are  constants 

The  first  strain  invariant,  I , is  defined  as  follows: 

e 

I = Zc  (17) 

c <x 

Geometrical  nonlinearity  is  introduced  into  this  relation  through 
the  strain,  which  we  choose  to  be  Seth's  n-measure  (23),  namely: 


cc  = -2_lL  (18) 

n 

Where  U }are  the  principal  stretch  ratios  n 

n is  an  empirical  constant  which  determines  the  measure  of 
the  strain 

An  additional  geometrical  nonlinearity  is  provided  by  the  dilatation 

term. 


A 


(19) 


Where  J (=  tiaJ  is  the  dilatation  ratio 

a . . 

6 is  an  empirical  constant 

Inserting  Equations  (17),  (18)  and  (19)  into  Equation  (16)  gives  the 
following  strain  energy  relation 


W = A): 


(20) 


Although  Equation  (20)  generates  a linear  stress-strain  law.  It  is 
both  materially  and  geometrically  nonlinear.  It  is  materially  nonlinear 
because  the  constants  A and  B (as  will  be  shown)  depend  Indirectly  on  the 
stress  through  the  porcsit/;  it  is  geometrically  nonlinear  because  of  the 
strain  and  dilatation  measures,  Equations  (18)  and  (IS). 
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Parameter  Determination 

The  parameters  A,  B,  n and  8 in  Equation  (20)  are  determined  by 
the  fact  that  W must  reduce  to  Hookean  behavior  in  the  limit  of  small 
strain.  We  proceed  as  follows: 

•.ln  = nAn  =n(l  + ne  ) = 1 + nl  * n2II  + n3III  (21) 

a a z c e 

£.  (i  +D 

J"6  = 1 - 8(1  + nil  +...)+  J!— — (n2I2+.. .)+... (22) 

e e 2 e 

When  Equations  (21)  and  (22)  are  introduced  into  Equation  (20),  we  obtain: 


W = AI 

e 

- BI  - Bn 1 1 + | ( 8+n) I2  + ... 
e e 2 e 

(23) 

In  the  limit  of  small 

strain,  we  have: 

0 = A - B 

(24) 

-2G  = -nB 

(25) 

K + = B(8  + n) 

(26) 

or 

A = 8 = r 

(27) 

2 

K-t  G 

3 

8 n 2G 

(28) 

WhereG  is  the  elastic  shear  modulus  of  the  propellant. 
K is  the  elastic  bulk  modulus  of  the  propellant. 
Thus,  our  strain  energy  function  is  gven  by: 
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Thus,  the  elastic  potential  contains  three  parameters  {n,  6,  and  K}. 
Such  a strain  energy  function  will  describe  the  elastic  behavior  of  many 
rubber-like  materials  reasonably  well  up  to  50%  strain,  or  thereabouts. 

The  actual  form  that  we  have  chosen  Is  unimportant  because  of  the  flex- 
ibility provided  by  the  parameter  n that  measures  the  strain  hardening 
or  strain  softening  response  of  the  rubber.  The  constant  n may  be  positive 
or  negative  and  Is  not  limited  to  integral  values. 

Now  the  fact  of  the  matter  Is,  solid  propellant  materials  evince 
hysteresis  upon  unloading.  One  way  of  explaining  this  is  that  the  two 
moduli  G and  K soften  during  loading  and  remain  fixed  during  unloading. 

In  order  to  represent  this  mathematically,  we  postulate  that  G and  K 
depend  upon  a parameter  y which  is  physically  related  to  the  area  fraction 
of  dewetted  particles,  and  which  is  controlled  by  the  Griffith  criterion. 
Our  picture  Is  as  follows.  In  a virgin  propellant,  we  have  zero  unboridlng 
or  zero  dewetting  as  shown  below: 
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no 

porosity 


After  loading , some  of  the  bonds  tear  and  dewet  and  we  have: 


load 


rubber  matrix 
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After  unloading  we  have: 


/ 


\ 


\ 


\ 


still 


no 

load 


porosity, 

unopened 


/ 


rubber  matrix 


¥ 


I 


t 


We  construct  a simple  model  as  follows.  Assume  the  unit  cell  to  be 
a sphere  of  radius  b containing  rubber  matrix  inside  of  which  is  a rigid 
filler  particle  of  radius  a.  The  contact  area,  <(> between  the  spherical 
particle  and  the  rubber  matrix  (whether  bonded  or  unbonded)  is  given  by: 


, 4ira 
* " 47.3 
3 b 
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where  $ Is  tl.e  area  per  unit  cell  volume. 

The  volume  fraction  of  filler  particles,  is  given  by: 


so  that: 


= K2/3 


(31) 


(32) 


Let  us  now  denote  the  fraction  of  the  area  that  Is  unbonded  by  <$ . 
Then  the  unbonded  area  per  unit  cell  volume  is  given  by: 


6<j> 


3 2/36 
bvf  6 


(33) 


If  we  denote  the  specific  surface  energy  associated  with  the  unbonded 
area  by  r,  then  the  total  surface  energy  of  the  unbonded  »egion  per  unit  cell 
volume  is  giv^n  by: 

% . 3 v2/3r6  (34) 

V+  b f 

Where  U3  is  the  total  surface  energy  of  the  sample 
V+  is  the  volume  of  the  undeformed  sample 

The  parameter  y (which  determines  the  softening  In  G and  K)  is  a 
measure  of  the  porosity  that  develops  when  an  unbonded  area  under  load 
opens  up  Into  a vacuole  as  shown  In  the  second  sketch  above.  Workers 
in  the  field  of  geophysics  (42)  have  shown  that  effective  moduli  can  be 
defined  in  terms  of  y in  the  form: 


and 


1 - Y 


Y 


(35) 


(36) 
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Equations  (35)  and  (36)  are  easily  derived  by  solving  the  elastic  problem 
of  a porous  sphere  under  load,  averaging  the  stress  field,  and  then  defining 
effective  moduli. 


Having  introduced  effective  moduli,  we  rewrite  Equations  ( 2$  and 
(29)  as: 


(37) 


and 


6 = n 


K - f Gn 
e 3 e 

2G_ 


We  now  postulate  that  y and  6 are  related  by: 

v = k6 


(38) 


(39) 


Where  k and  m are  constants 

Equation  (39)  states  that  the  porosity  y which  occurs  under  load  is  a 
nonlinear  function  of  the  area  fraction  of  unbonded  matrix  6.  Because  the 
actual  determination  of  the  porosity  via  a stress  analysis  is  an  extremely 
difficult  problem  and  because  of  interactions  between  adjacent  cells,  we 
prefer  to  eschew  this  difficult  analysis,  and  introduce  one  extra  parameter 
m,  an  interaction  parameter,  into  the  theory.  The  proportionality  constant 
k will  turn  out  to  be  absorbed  into  the  fourth  parameter,  which  is  the  surface 
energy. 


Stress-Strain  Relations 
The  stress-strain  law  is  given  by: 


_ oW_  _ __e  / n _ ,-B\ 

A a = o J = A oA  n 'a  ' 

u a a a a 


(40) 


where 


oa  is  the  principal  Piola  stress 


and 


oq  Is  the  principal  Cauchy  stress 
The  power  input  to  a sample  is  given  by: 


a A 
a a 


(41) 
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The  total  strain  energy,  , is  given  by: 

U1  = V+W  (42) 

The  total  surface  energy  is  given  by: 

u3  = v+  l J*  r(*)1/m  (43) 

and  t»ie  Griffith  criterion  postulates  that: 

+ U3  (44) 

where  the  dots  denote  time  derivatives. 

We  now  turn  to  the  case  of  simple  tension  with  superimposed 
hydrostatic  pressure.  For  this  case,  we  have: 


X-J  = X 

1 /IC\ 

X2  = X3  = JT 

(46) 

o ] = S - P 

(47) 

o2  = o3  = -P 

(48) 

Upon  substituting  into  Equation  (40),  we  obtain: 


2G„ 


SO 


$ 


n/2 


(49) 


2G 


PJ 


e 

n L 


O'*-  (i) 


n/2 


(50) 


After  some  simplification,  Equation  (44)  reduces  to: 


40 


n(4+y)‘ 


l n/2 

>n  + 2(t)  -3  J'g-  1 

n + 3 


16 


36y2  L 


2/3 


1-m 


3r  y m 


(51) 


bG  mk 


1/m 


Equations  (49),  (50),  and  (51)  complete  define  S and  P In  terns  of  x, 
J,  and  y. 
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EXPERIMENTAL  EVALUATION  OF  CHARACTERIZATION  PARAMETERS 
For  this  effort  we  will  define  a new  constant  M as  follows. 


M = 


rv 


2/3 


bG  mk 


1/m 


(52) 


Thus,  the  basic  material  characterization  requires  the  evaluation 
of  four  parameters  n,  G,  m and  M.  The  problem  is  to  work  backwards  from 
available  data  and  determine  these  four  parameters  in  a self-consistent 
fashion.  First  we  set: 

A = ee  or  e = £nx  (53) 

J = eA  or  A = £nJ  (54) 

whereupon  Equation  (49)  may  be  rewritten: 


1 


n S a 
2 G e 


1 +4  G 


,ne 


- e 


-n 


e-A 

2 


(55) 


Likewise  Equation  (50)  may  be  ^written: 

r 
en 

4 - v 


n 


6 y 


= 6 = 


n *2  + | (e"E 


- e'n  T") 


e-A 


(56) 


We  now  consider  two  cases.  First,  take  the  case  P - 0.  Since  the 
strains  we  are  dealing  with  do  not  exceed  50%,  we  can  expand  Equations  (55) 
and  (56)  to  second  order  and  obtain: 


1 + 


ne 


A_ 

3e 


(57) 


and 


(58) 
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When  we  eliminate  y between  Equations  (57)  and  (58),  we  obtain: 

S _ En  , r / e - 2a  \ 

E[~r) 


(59) 


rt  Pa 

Equation  (59)  suggests  that  a plot  of  ($/e  ) vs  should  yield  a straight 

line  with  slope  E and  Intercept  E n/4.  In  addition,  when  we  expand  Equation  (51  ) 
we  obtain: 

1 m 1+m 

1-m  "1+m  1-m  (60) 

M e 


*■  (f) 


which  suggests  that  a log-log  plot  of  4 vs  e should  determine  M and  m via 
slope  and  intercept. 

Atmospheric  Pressure  Tests 

Figures  46  to  50  display  the  ] og-1 og  correlations  of  a vs  e for  Farris' 
test  data  #1,  5,  6,  7,  8,  and  10  at  e = 0.5  min"  , #2  at  0.25  min-';  §3  at 
0.05  min-1;  #4  at  0.0125  min"l;  and  #9  at  1.25  min^.  In  all  cases,  good 
straight  lines  are  obtained.  The  slight  departure  from  linearity  at  the 
higher  strain  levels  arises  from  the  fact  that  we  approximated  the  fracture 
criterion  by  way  of  second-order  expansions.  If  we  had  kept  the  complete 
function  and  put  it  on  the  computer,  we  could  have  developed  complete  linearity. 
For  the  test  cases  enumerated  above,  we  find: 


TABLE 

II 

Parameter  Evaluations 

from  Tests  at  77°F 

and  0 

psig 

e,  min"1 

0.0125 

0.05 

0.25 

0.59 

1.25 

m 

0,55 

0.54 

0.51 

0.56 

0.57 

M 

0.070 

0.079 

0.089 

0.091 

0.076 

<m>  = 0.55 


Thus  the  fracture  parameters  m and  M are  remarkably  constant  and  indepen- 
dent of  strain  rate. 


g-2a 


Figures  51  to  55  display  the  linear  correlation  between  (S/e^)  and 
for  the  same  five  data  sets.  For  these  cases,  we  find: 
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P = 0 psig 
e = 0.5  min.  * 


#1 

#5 

#6 

#7 

#8 

#10 


O Test  # 4 


Figure  50,  Fracture-  controlled  Dilatation  of 


Aerojet  CTPB  Propellant 


Figure  52,  Fracture-  controlled  Simple  Tensile  Behavior  of  Aerojet  CTPB  Propellant 


Fracture- controlled  Simple  Tensile  Behavior 
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TABLE  III 


Parameter  Evaluations  from  Tests  at  77°F  and  0 psig 


e,  min"^ 

0.0125 

0.05 

0.25 

0.50 

1.25 

E,  psi 

385 

400 

567 

555 

692 

n 

1.04 

2.00 

1.41 

1.44 

1.16 

Note  that  E increases  monotonlcally  with  increasing  strain  rate.  This 
trend  represents  the  effect  of  viscoelasticity  which  we  have  left  out  of  the 
theo\\  One  should  not  read  too  much  significance  into  che  scatter  of  n 
because  n is  determined  from  the  intercept  of  a line  in  such  a way  that  a 
slight  change  of  slope  produces  a large  change  in  n.  This  sensitivity  arises 
because  the  plot  involves  the  large  numbers  (S/e2).  This  scatter  will  be 
markedly  reduced  when  the  whole  program  is  put  on  the  computer.  In  this 
initial  cut  at  the  theory,  we  reduced  all  the  circulations  to  an  approximate 
form  which  enabled  us  to  make  log-log  or  linear  plots. 

Superimposed  Hydrostatic  Pressure  Tests 

Now  we  consider  the  case  where  P>  0.  Looking  at  Equation  (56)  we  note 
that,  when  P >3S  the  numerator  of  Equation  (56)  becomes  negative,  while  A 
remains  zero  or  positive.  This  reflects  the  bi -nonlinear  aspects  of  propellant 
behavior,  fint  clearly  set  forth  by  Dong,  Herrmann,  Pister,  and  Taylor 
(43).  Putting  It  another  way,  if  the  filler  particles  were  not  present,  the 
rubber  would  contract  under  pressure  and  A would  be  negative.  But,  because 
both  the  rubber  and  Filler  particles  are  incompressible  the  hydrostatic  behavior 
of  propellant  is  discontinuous,  and  a plot  of  a vs  P looks  like  the  following 
sketch. 


Thus  Equation  (56)  must  be  replaced  by  the  statement  a = 0,  when  P > S/3. 
When  we  do  this,  we  find  the  following  relation: 


S 

T 

e 


1 - 


1 

r 1— m 

<f>  <y 


2 


m 

1-m 


= E 


- + E( 
4 


c - A/3\ 
2 1 


(61) 
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13  U WKt0?!nal  Vhe  Hu  au  50  pstg  and  °*5  min'1  strain  rate  (Test  #'s 
*5»  16)  and  used  the  M and  m determined  from  fho  *1  n • * s 

and  made  the  straight  line  correlation  suggested  by  (61)  Thp  Lc  ?+Tf 

shown  in  Figure  56.  The  agreenent  is  fantl!tic  We  get  ihe 

lUJe  nS=1!4™’ned  a*  “r0  PreSSUre-  We  9et  *•»  E - 555.  |Si! 

100  u"3S5  r57t0  ZsTi  £ 

line  is  excellent,  and  the  modulus  (E  = 550  ns  H*  fn*  *9?  " strai9ht 

fin  :■  sr  ? •P7 

thlcTl s hextfrere ly ' lensl t? Je If  the ’s “pe  f 'thfft "1  P"’""*"*’* 


50  psig 
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MODIFICATION  OF  THE  THEORY 

As  shown  above  the  previous  theory  provided  excellent  correlations 
with  uniaxial  test  data  for  pressures  up  to  50  pslg.  The  purpose  of  this 
effort  was  to  extend  the  correlation  up  to  500  pslg. 


In  1966,  Tobolsky  and  Shen  (44)  suggested  that  the  strain  energy 
function  for  NR  be  modified  by  the  factor  J*  to  explain  the  dependence  of 
the  force-pressure  coefficient  upon  strain.  Following  this  Idea,  which  has 
been  used  successfully  by  Sharda,  we  write  the  strain  energy  factor  In  the 
form: 


2Ga  « *n-l  i-e  ^ 


(62) 


where  « is  a constant 

For  pressures  above  600  pslg  the  dilatation  Is  essentially  suppressed; 
which  gives: 


Ja  = 1 , P > 600  psig  (63) 

Under  these  conditions  Equation  (62)  reduces  to  Equation  (37).  Hence 
the  parameter  characterization  results  given  In  Tables  II  and  III  remain  unchanged. 


Following  the  procedure  outlined  In  the  preceding  section,  and  omitting 
unnecessary  algebraic  details,  we  arrive  at  the  following  expression  for  simple 
tension  under  superimposed  hydrostatic  pressure,  after  expanding  to  second 
order: 


Y 


(64) 


- Is-  <r * '~sr  - TO 2 Ti  + 3«>] 


(65) 


Figure  59  shows  a plot  of  the  zero  pressure  data,  from  which  we  find: 

m = 0.60  (66) 

M = .074  (67) 


£+  3a  = 2.0 


(68) 


The  value  of  n previously  determined  from  the  200  psi  data  was  5.33. 
Therefore  we  calculate  that: 

a = -0.22  (69) 

which  is  In  remarkably  good  agreement  with  the  value  determined  by  Sharda  (45,  46) 
and  by  Tobolsky  and  Shen  (44),  namely  p = -0.20. 
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Introduction 


In  this  report  we  present  a rationalization  of  the  observe!  fact  that, 
in  propellants , stress-at-break  seems  to  correlate  with  the  squire  root  of  the 
relaxation  modulus  evaluated  at  time-to-break.  The  basis  for  this  rationaliza- 
tion is  the  ablating-spherical-hole  model  of  Max  Williams,  which  has  never 
been  correctly  analyzed  for  reasons  which  we  shall  present  below.  The  beauty 
of  the  model  however  lies  in  the  observation  that  the  spherical-flaw  model  gives 
basically  the  same  fracture  prediction  as  does  the  penny -shaped-flaw  model, 
the  two  differing  only  in  a shape  factor.  The  advantage  of  the  model  is  that  the 
algebra  associated  with  its  development  is  relatively  simple.  The  basis  for 
using  it,  or  the  justification,  is  obtained  if  one  imagines  that  the  applied  stress 
(to  a simple  tensile  bar)  is  concentrated  around  a void  in  some  fashion  that 
makes  the  void  grow  like  a bubble  expanding.  This  stress  concentration  effect 
means  that  our  final  results  on  the  spherical-flaw  must  be  related  back  to  the 
observed  data  on  a tensile  bar  simply  by  multiplying  in  a stress  concentration 
factor.  This  is  actually  a somewhat  crude  approximation,  which  will  be  late.* 
replaced  (and  discussed  in  a subsequent  progress  report)  by  a more  sophistica- 
ted randomly-flawed  model  which  is  currently  under  development. 

Analysis 

In  order  to  set  the  theme  fcr  this  discussion,  we  must  understand  the 
general  theory  of  viscoelastic  constituti  ve  equations . In  this  repost  we  restrict 
the  discussion  to  linear  behavior.  The  linear  theory  of  viscoelasticity  is  con- 
veniently represented  by  a Max/zell  model,  which  is  a parallel  array  of  one 
rubbery  spring  and  a distributii. -n  of  Maxwell  elements;  each  element  is  a series 
array  of  one  Hookean  spring  ar.d  one  Ne  wtonian  dashpot;  the  Maxwell  elements 
nre  distributed  in  terms  of  their  relaxation  times.  The  rules  for  the  model  are: 
stress  adds  in  parallel  and  is  constant  in  series;  strain  adds  in  series,  and  is 
constant  in  parallel.  With  this  in  mind,  we  have  the  following  relations: 


t * tR+S  t1 

(A-l) 

tR-2GR.t+  (*R-§  Gr)  * l 

(A-2) 

(A-3) 

i . ID  , » . ID 

i + (KiVf  Vi)”  i 

(A-4) 

( A-5) 
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where 


t is  the  total  stress  tensor 


is  the  stress  tensor  in  the  rubbery  spring 
t*  is  the  stress  tensor  in  any  Maxwell  element 
e is  the  applied  strain  tensor 


I is  the  idem -tens or 


G and  K are  the  elastic  moduli 


are  the  shear  relaxation  times 


fB1  are  the  bulk  relaxation  times 

IS 

e is  the  3train  tensor  in  the  spring  of  a Maxwell  element 


e*^  is  the  strain  tensor  in  the  dashpot  of  a Maxwell  element 
What  we  do  now  is  to  separate  the  above  expressions  into  traces  and  de viators 
We  have: 

t R » sR  - p8! 


■**+ 


i i T 
* S - p I 


(A-6) 

(A-7) 


e * ? + $ I 

~ ~ T*  ~ 


(A-8) 


,HS  , ,1S  , 9iS  T 

•3  * C + v I 

~ - -3 — 


(A-9) 


e10  * <iD  + «?1D  I 


(A-10) 
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where 

s is  the  stress  deviator 
c is  the  strain  deviator 


From  (A-2),  we  have: 


& 

1 

K 

(A-n) 

sR  * 2GC  < 
~ R ~ 

(A-12) 

From  (A-3),  we  have: 

p*  * _Kjl 

(A-13) 

s1  *=  2Gi  € iS 

(A-14) 

From  (A-4),  we  have: 

* - -K  T ^ 

P *i  TBi 

(A-1S) 

s1  = 2G  ^ 

(A-16) 

From  (A-5) , we  have: 

*iS  + - * 

(A-17) 

f1S  + €1D  * € 

(A-13) 

A>  /V 


When  Equations  (A-13),  (A-lb)  and  (A-17)  are  combined,  we  obtain* 

* lo  e " dW  ^A_19) 


(A-20) 

TB1 

i 

P = 

-K,  tf1S 

(A-21) 

Likewise,  when  Equations  (A-14),  ( A- 1 6 ) , and  ( A- 1 8 ) are  combined,  we  obtain: 
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IS 


t^ 

ri. 


J*  e~  t4  d € (t?) 


(A-22) 


,iS 


i® 


s1  * 2G,  €iS 

A/  i <V 


When  we  add  everything  up,  we  have: 


(A-23) 

(A-24) 


♦ t - 11  t 

l = 2f  (C^+SGi  e ~ — ) d c (f)  +2  J (KR+SKie"rB1)d  * {t') 
o o 


( A-2E ) 


In  order  to  simplify  the  ensuing  analysis  in  a nonessential  fashion  , we 


shall  assume  that  all  of  the  relaxation  processes  are  shear  processes,  and  that 
the  bulk  processes  are  purely  elastic.  This  implies  that  the  7U.  are  all  zero. 
Then  Equations  (A-19,  A-20,  A-21  and  A-25)  become:  1 

(A-26) 

i'K  o 

(A-27) 

P1  = - K,  * 

(A-28) 

t-t  ' 

_t=2f  (G^EGj  e 
o 

(A-29) 

where 

Kg  - 

(A-3C) 

The  stored  energy  in  the  model  is  the  sum  of  the  strain  energies  in 
all  springs.  In  any  one  spring,  the  strain  energy  is  given  by: 

Wl=4ti:e1S 

4 A#  A» 

(A-31) 

-iV-P‘i).<«“  ,£n 

4 Vv  ~ ~ 

(s1:  <1S-  p1  ^1S) 

**  A/  r 

(A-32) 

.G,  <iS.ciS  + Kl  (*IS>2 
i ~ ~ "o' 

(A-33) 
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The  stored  energy  is  thus: 

W ■ O,  c : € + S G,  €aS  : €iS  + KG  *2 

K ~ ~ i~  ~ ~2 — 

likewise,  the  dissipation  rate  is  given  by: 


Et1:  S10  = E (s^p1 


, , *1S  s1:  P 


2E^i-<1S;  <lS 

<r  <v 

Ti 


The  spherical  flaw  model  deals  with  an  ablating  hole  of  radius  a in 
a sphere  of  radius  b,  subjected  to  isotropic  tension  at  r = b.  The  displacement 
must  have  the  form: 


u « Ar  4 


whereupon 


eQ  A + "]X 

. a , 2B 

er°A'T 

r 


* 3A 

€ = JL 

ce  r3 

c = - 2B 
r r3 


o = - 


4 pt 
3 J 

r o 


J*  (Gr  + SGje  ~ Tj  ) dB(t')  + 3KqA 


00*  “tJo  (%i'SGie  Ta  ) dB(t')  + 3KqA 


The  radial  stress  is  zero  at  r = a,  giving: 


so  that: 


3KGA  = "X  J*  ^GR+  SGie  ri  ) dB^') 


a o 


B +ir  / (G  tSG  e-  V ) dB(f) 
r 3K.-  az  o 1 
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r “ 4(“3  " ~3>  J*  (Gfc+  SG±  e ’ Tt  ) d B(f) 
a r o 

°8  !S(f3+^)^o  «VSV"^  dBft‘) 


At  r - b#  the  radial  stress  is  S,  giving: 


S * 4(J5  " ~3}  Sq  (Gr+  ^e  “ ) dB(t‘) 


The  power  input  to  the  sphere  is  given  by: 


S • 4yb2  • u, 


where 


•*  3 4b  a <»*  t-t* 

ub  = TF  “ 4 J (Gr+  SG  e ' r,  ) dB(t') 

b Kq  a o K 1 1 

+.4b r * 4b  r*  Gi  - — 

2K  a3  GG  B ~ r 3l  dD<1)  2T~  e T 

3KGa  3KGa  o Ti  1 


Let  us  define: 
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*i“J  e Ti  dB(t') 
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The  stored  e nergy  is  given  by: 


W * (GrB2  + L Gi<pi  2)  + g-  (grb  + sg1 


KGa 


The  total  strain  energy  in  the  sphere  is  given  by: 


„ r>  2.  Wb3 

Uj-  4 W4irrdr»— (-j--  -jj  ^ 


(A-55) 


(A-56) 


The  dissipation  rate  is  given  by: 


®=1f 

r6  T1 


(A-57) 


And  the  total  dissipation  rate  is  given  by: 


U»  * 16 


TT  A JL\  £ Gl-S- 

V b3J  Ti 


The  total  surface  energy  is  given  by: 

U3  « T 4na2 

And  the  Griffith  criterion  demands  this: 
9'  Uj  + U2  + U3 


which , after  simplification , leads  to: 
0 


* r*5 


» c-^ ( GrB2  + £ G^2) 


^g(GrB*  EG^/j 


(A-58) 


(A-59) 


(A-60) 


(A-61 ) 


The  third  term  In  the  brackets  of  (A-61)  can  be  neglected,  whereupon  (A-61) 
becomes : 

5 

0 - a [^f (GrB2  + S Gi^»i2)  ] (A-62) 

Equations  (A-53)  and  (A-62)  jointly  relate  (a,  B,  and  Z) . Equation  (A-62)  explicitly 

--  2 2 

states  that  the  void  does  not  grow  (a  * 0)  until  [ — = GRB  + JIG^  ]. 

This  defines  the  criticality  condition  at  which  fracture  starts. 
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We  are  interested  in  relating  S to  a and  also  determining  the  time-to- 
break,  when  ( a »b)l  B can  be  eliminated  from  (A-53)  by  Laplace  transformation. 
However,  there  is  an  easier  way  to  carry  out  the  analysis  without  getting  in- 
volved in  the  solution  of  a nonlinear  integral  equation.  We  assume  that  B in- 
creases exponentially  with  terms  of  the  form: 


B **  Ae 


at 


Then  (A-53)  becomes: 


a'  b 

In  writing  (A-64),  we  have  assumed  that: 
eat>>  e'^l 


Tr'rn  -b(gr+egi  t^) 


(A-63) 

(A-64) 


(A-65) 


There  are  two  cases: 

a)  if  t is  greater  than  , then  (A-65)  is  good  for  any  a 


b)  if  t is  smaller  than  then  (65)  is  good  for  at  > 1 

When  (A-63)  is  substituted  Into  (A-62),  we  obtain: 

r _5  9 Ct  T.  ? 

-T-  ’B ‘cg**  igi(tt^-)!] 

Elimination  of  B yields: 


S.4(l--si)(GR+SGt-^r-) 


T/3a 


(A-66) 


(A-67) 


GR  + EG 


pL— \ 

i\l  + a tJ 


Equation  (A-67)  determines  A as  a function  of  a,  or  vice  versa.  For  the  purely 
elastic  fracture,  we  have: 


S = 4 


(A-68) 


so  that  (A-61)  really  determines  an  effective  modulus: 
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The  parameter  a is  determined  by  the  break-time  t^  which  in  turn,  according 
to  the  model,  is  determined  when  a = b.  Equation  (A-66)  yields: 


if  = A2  e2* 


*b  GR+EGi(r^)] 


Let  us  now  consider  each  of  the  two  limiting  cases: 

a)  when  t > we  take  a < < 1,  and  arrive  at: 


where 


where 


(l^^) 


eff  R 8b  \3A2Gr  y 


ti  * L Gt  Ti 


b)  when  t < Timin,  we  take  a > > 1 , and  arrive  at: 


-ai  .A2e  2atb  Gg  or  ax  t *(b£- 

3 G 2tb  Wg, 


Ge££  " GG  + 


r y v 2 i 

i ( GiV  Gi 

~2~  ~z~  ) ~ Gp  E — y 

„ a2  \ V G r2 

a «-  1 J 


GG  “ gr  + L Gi 


(A— 70) 


(A-71) 


(A-72) 


(A-73) 


(A-74) 


(A-75) 


(A-76) 


(A-77) 


Por  intermediate  values  of  a or  t,  , (A-69)  and  (A-70)  can  be  put  on  the  computer 
for  an  assumed  distribution  function  G4 . A very  useful  representation  is: 


A-9 
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E 


1 4-aTj 


1 

1 + k an 


(A-71) 


where  0 < n < 1.  This  is  the  Cole  distribution  function. 

In  summary,  what  we  have  shown  is  that  the  stress-at-break  is 
proportional  to  the  square-root  of  the  effective  relaxation  modulus-at-break 
(A-68) , where  the  effective  modulus  Is  defined  by  A-69),  and  the  tlme-to-break 
Is  Introduced  through  (A-70). 
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EFFECTS  OF  PROPELLANT  DEWETTING 


The  oxidizer/binder  dewetting  phenomenon  has  a number  of  effects  on 
propellant  response  and  failure  behavior.  The  following  sub-sections 
describe  some  of  the  major  effects;  namely:  (1)  The  effect  of  dewetting 

on  propellant  failure  criteria;  (2)  the  effect  of  superimposed  pressures 
upon  vacuole  formation;  and  (3)  the  dependence  of  low-strain  moduli  upon 
dewettinq. 

1.  The  Effect  of  Dewettinq  on  Propellant  Failure  Criteria 

Previous  studies  (B-l  to  B-5)  at  Aerojet  produced  a simple  fail- 
ure criterion  for  solid  propellants.  This  criterion  was  developed  from 
constant  load-to-failure  testinq  usinq  uniaxial,  strip-biaxial  and  poker 
chip  tests.  In  its  simplest  form,  all  test  data,  regardless  of  testino 
mode,  can  be  reduced  to  a single  curve.  This  curve  represents  the  time- 
to- failure  under  a constant  applied  stress.  A typical  failure  curve  is 
shown  in  Figure  B-l.  The  curve  was  obtained  using  data  from  two  different 
cartons  of  a polybutadiene  propellant  (using  the  ay  values  obtained  by 
Lockheed  Propulsion  Company  on  their  STV  program).  As  shown  in  Reference  B-l, 
test  data  from  strip-biaxial  and  poker  chip  tests  fall  on  this  failure  curve. 

The  same  curve  was  obtained  for  tests  performed  under  superimposed 
hydrostatic  pressures  (Fiaure  B-2).  However,  in  this  case  it  was  necessary  to 
shift  the  test  data  along  the  time  axis  by  a time-pressure  shift  factor,  ap. 

A plot  of  ap  versus  pressure  is  provided  in  Figure  B-3.  These  results  are 
typical  of  those  obtained  for  many  other  propellant  formulations. 


B-l 


EQUIVALENT  CONSTANT  STRESS  TO  FAILURE  CURVE  OBTAINED  FROM 
TWO  CARTONS  OF  THE  STV  PROPELLANT 


min, 


Figure  B-2.  EFFECT  OF  SUPERIMPOSED  PRESSURE  ON  CONSTANT  STRESS  TO  FAILURE  DATA 

FOR  THE  STV  PROPELLANT  (PRESSURE  SHIFTED  TO  0 PSIG) 
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The  siqnificance  of  these  results  is  best  seen  when  they  are 
plotted  in  principal  stress  space.  The  failure  surface  so  generated  is 
not  that  predicted  by  any  of  the  classical  theories,  yet  it  is  straight- 
forward and  easily  interpreted.  As  a convenience  in  presenting  the  data, 
we  shall  consider  the  results  in  two-dimensional  stress-space.  The  curves 
presented  in  Figure  B-4  are  crossplots  from  Figures  B-2  and  B-3  taken  at 
discrete  levels  of  time- to- failure  under  constant  stress  (from  0.01  min.  to 
1 year).  Each  curve  represents  the  propellant  failure  criterion  for  a 
given  time-to-failure. 

The  data  in  the  first  quadrant  (where  both  principal  stresses 
are  tensile)  are  typical  of  the  classical,  maximum  principal  stress  failure 
criterion.  It  is  considered  reasonable  that  this  should  be  the  case,  since 
the  propellant  is  hiqhly  dewetted  prior  to  failure.  Hence,  the  lateral 
stresses  should  have  only  minor  effects  upon  the  behavior  of  the  propellant. 

Ir  -«e  second  and  third  quadrants,  a surprising  effect  is  observed. 
Under  small  compressive  stresses,  the  strength  of  the  propellant  actually 
increases.  This  is  clearly  evident  in  the  curve  for  the  time-to-failure 
equal  to  0.01  min.,  where  it  is  seen  as  a marked  "shoulder"  on  the  failure 
surface.  This  shoulder  occurs  at  each  of  the  failure  times,  but  is  less 
evident  at  the  lonaer  times.  The  shoulder,  in  Figure  B-4,  is  believed  to 
occur  because  of  an  interaction  between  the  dewettina  process  and  the 
superimposed  hydrostatic  pressure.  In  the  tensile  quadrant,  the  propellant 
is  extensively  dewetted  before  failure.  In  the  second  quadrant,  the  addition 
of  hydrostatic  pressures  acts  to  suppress  dewettina.  Extension  of  this 
arqument  suggests  that  dewetting  will  be  suppressed  with  increasina  hydro- 
static pressure  until,  at  least  asympotically,  the  dewetting  will  be  com- 
pletely eliminated  prior  to  propellant  failure. 
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By  this  sort  of  reasoning,  it  is  possible  to  define  the  maximum 
uui axial  strength  that  the  propellant  might  have  achieved  if  dewetting 
had  not  occurred.  As  shown  in  Figure  B-5,  extrapolation  of  the  estimated 
pressure-defined  asymptote  back  to  0 psig  pressure  uives  a uniaxial 
strength  of  340  psi  for  the  undewetted  propellant.  The  dewetted  propellant 
was  found  to  have  a strength  of  only  224  psi,  which  is  about  66%  of  its 
potential  strength.  Thus,  dewetting  decreased  the  strength  of  this  propellant 
about  34%  for  a time  to  failure  of  0.01  min.  The  percentage  reduction  in 
strength  was  found  to  be  about  the  same  at  all  times- to- failure,  see  Table  B-l. 

In  this  section  we  have  presented  a propellant  failure  criterion 
with  one  possible  interpretation  for  the  mechanisms  involved.  The  next 
section  presents  some  data  which  sugaest  a different  mechanism  for  the 
processes  involved. 

2 . The  Effect  of  Superimposed  Pressures  on  Vacuole  Formation 

Farris  (B--6,B-7)  has  shown  that  the  effect  of  superimposed  pressures 
is  to  delay  the  onset  of  vacuole  formation.  Figure  B-6  was  reproduced  from 
Reference  3-6,  to  demonstrate  how  increasing  pressures  appear  to  suppress 
the  dewettinq  processes.  An  anomaly  exists  here,  since  the  suppression  of 
dewetting  would  be  expected  to  make  the  propellant  behave  more  rigidly  and 
exhibit  marked  reductions  in  failure  elongation.  However,  the  opposite 
seems  to  be  true,  the  elonqation  at  break  increasing  with  increasing  pressure. 
Thus,  the  failure  process  appears  to  be  one  where  (a)  dewetting  occurs,  but 
because  of  the  superimposed  pressure,  the  vacuoles  are  suppressed  and  (b)  the 
stretching  mechanism  of  the  propellant  is  restricted  to  local  shearino  Dro- 
cesses. 
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TABLE  B-l 

COMPARISON  OF  OBSERVED  PROPELLANT  DEWETTED  STRENGTH  WITH 
THOSE  ESTIMATED  FOP  THE  UNDEWETTED  MATERIAL 


Observed  Dewetted  Strength 


Time-to-Failure 

Estimate  of 
Undewetted  Strength, 
psi 

Strength, 

psi 

Percent  of 
Estimated 
Strength 

.01  Min. 

340 

224 

66 

1 Min. 

192 

130 

68 

1000  Min. 

86 

58 

68 

1 Year 

41 

27.6 

67 
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COMPARISON  OF  THE  STRESS-STRAIN  AND  DILATATION-STRAIN 
BEHAVIOR  OF  A HIGHLY  FILLED  ELASTOMER  (65  Vol.  %)  AT  A 
SERIES  OF  HYDROSTATIC  PRESSURES 


In  summary,  it  appears  that  the  propellants  may  fail:  (1)  by  a 
maximum  principal  stress  failure  criterion  at  ambient  pressure;  and  (2) 
by  a shear  failure  criterion  under  hydrostatic  pressures.  The  transition 
from  one  criterion  to  the  other,  as  described  earlier,  may_ be  restricted 
to  a fairly  narrow  range  of  pressures. 
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